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ABSTRACT 

Observations of strong gravitational lensing, stellar kinematics, and larger-scale tracers enable ac¬ 
curate measures of the distribution of dark matter (DM) and baryons in massive early-type galax¬ 
ies (ETGs). While such techniques have been applied to galaxy-scale and cluster-scale lenses, the 
paucity of intermediate-mass systems with high-quality data has precluded a uniform analysis of 
mass-dependent trends. With the aim of bridging this gap, we present new observations and analyses 
of 10 group-scale lenses at (z) = 0.36 characterized by Einstein radii 0Ein = 2"5 — 5"1 and a mean 
halo mass of M200 = 10^"^ ° M©. We measure a mean concentration C200 = 5.0 ± 0.8 consistent with 
unmodified cold dark matter halos. By combining our data with other lens samples, we analyze the 
mass structure of ETGs in 10^^ Mq — 10^® Mq halos using homogeneous techniques. We show that 
the slope of the total density profile ytot within the effective radius depends on the stellar surface 
density, as demonstrated previously, but also on the halo mass. We analyze these trends using halo 
occupation models and resolved stellar kinematics with the goal of testing the universality of the DM 
profile. Whereas the central galaxies of clusters require a shallow inner DM density profile, group- 
scale lenses are consistent with a Navarro-Frenk-White profile or one that is slightly contracted. The 
largest uncertainties arise from the sample size and likely radial gradients in stellar populations. We 
conclude that the net effect of baryons on the DM distribution may not be universal, but more likely 
varies with halo mass due to underlying trends in star formation efficiency and assembly history. 
Subject headings: dark matter — galaxies: elliptical and lenticular, cD — gravitational lensing: strong 


1. INTRODUCTION 

The distributions of dark and baryonic matter within 
galaxies of various masses is a key constraint on theo¬ 
ries of galaxy formation and cosmology. In the stan¬ 
dard cold dark matter (CDM) model, the distribution of 


well understood from Wbodv simulations (e.g 

, Navarro 

et al. 

1996b| Diemand et al. 

2005 

Gao et al. 

2U12|). In 


onic physics, however, the distributions of stars, gas, 
and DM depend on poorly understood processes such 


20041), thermal and mechanical feedback from supernovae 

iJNavarro et al. |1996a 

Pontzen & Governatol 

2012 

) and 

active galactic nuclei (AGM; Martizzi et al. I2U13 

, and 

dynamical heating in mergers (e.g., El-Zant et af. 

2001 

Nipoti et ai.||2004l ITonini et ai.||20U6| (Laporte & 

White! 

2015|). 

Detailed observations of the mass distribution there- 


fore contain important information on the balance of 
these competing baryonic processes. Of particular in¬ 
terest is the radial density profile of DM on small scales, 
which is sensitive to this balance and may als o constrain 


the microphysics of the DM particle (e.g., Spergel & 
Steinhardt 20001. Since the relative importance of the 


various baryonic processes is expected to vary with a 
galaxy’s mass and formation history, a valuable route to 

^ Present address: European Southern Observatory (ESO), 
Karl-Schwarzschild-Strasse 2, 85748 Garching, Germany 


progress is to examine the distributions of dark and bary¬ 
onic matter across galaxies, groups, and clusters, thereby 
spanning the full range of systems where the relevant ob¬ 
servational techniques can be applied. 

Strong gravitational lensing has emerged as a key tech¬ 
nique for tracing the mass distribution for this wide range 
of systems, since it provides a geometric mea sure of the 


tota l mass within the E instein radius (see, e.g., Treu 2010|| 
and Treu fc Ellis||2014 for recent reviews). For the more 


massive systems, weak lensing allows the total mass to 
be traced to larger scales. Other observations, such as 
stellar kinematics on smaller scales where the stellar con¬ 
tribution is usually dominant, and satellite dynamics and 
X-ray emission on larger scales, enable the mass profile 
to be measured at several widely separated radii. This is 
essential to constrain multi-component models that sep¬ 
arate the stellar and DM components. 

The combination of strong lensing and stellar kinemat- 


of early-type galaxies (ETGs; 

Treu & Koopmans 

2002 

2004||Jiang & KochaneK|2007l 

Auger et al.|2010al| 

van de 

Ven et al.j 20101 jBarnabe et a 

1. 20111 2013 Grlllo et al. 

2013j). Based on more than 100 lenses discovered in the 


the logarithmic slope 7tot (also denoted 7') of the to- 
tal density profile within the effective radius i?e, where 
Ptot oc has a mean value (7t,nt ) = 2.078T0.027 an d 

a fairly small scatter of 0.14 ± 0.02 (Auger et al.||2010a). 
This has been interpreted as evidence tor a “conspiracy” 
between DM and baryons that drives their combination 
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to a nearly isothermal density profile (Koopmans et al. 

2006 

2009 

Treu et al. 

2006 

Gavazzi et al. 2007). 'i’he 

SLAC 

S sample and more recent surveys (SL 2 S: Kutt et al. 

2011 

Sonnenfeld et al. 2013a 

BELLS: 

Bolton et al. 2012|) 


components, but to date o nly two group-scale lenses have 


tion (IMF) of massive ETGs (Auger et al. 2010b Treu 


et al. 12010 ) as well as the evolution ot the density profile 


over cosmic time (Sonnenfeld et alT]|2013b 2014). All of 


these studies pertain to galaxy-scale lenses with Einstein 
radii in the range 1 " < 0 Ein ^ 2 ", velocity dispersions 
cr « 250 ±40 km s~^, and halo masses M 200 — 10^^’^ Mq 
( G avazzi et al.|[2007 |. 


Similar techniques have been extended to giant ellip¬ 
ticals in the centers of massive, relaxed clusters with 


iM 200 - 
2002 

2004 

2008 

) and 


which cu lminated in a study ot 7 systems 


by sand et al . 
i2009i I2011D, 

Newman et al. 


2013a|b hereafter N13a, N13b, or N13 collectively], the 


average total density slope within i?e was found to be 
(7tot) = 1.16±0.05]']g Q7. This is much shallower than for 
galaxy-scale lenses and consistent with high-resolution 
dark matter only simulations, despite the presence of sig¬ 
nificant stellar material on these scales. As a result, after 
decomposing the density profile into its stellar and DM 
components, N13 found the inner DM density profile to 
be shallower than the canonical NFW slope, an intrigu¬ 
ing result which several group of simulators have sought 
to explain ( Mar tizzi et al.||2013 Laporte fc Whit^|20i5 


Schaller et aT 


. |2015bt . 
il questioi 


i'wo natural questions are the origin of the wide range 
of total density profiles seen in ETGs, and whether the 
shallow DM profiles evident in brightest cluster galax¬ 
ies (BGGs) are also present in lower-mass systems. The 
main goal of the present paper is to connect the trends 
observed in galaxy- and cluster-scale lenses by consider¬ 
ing ETGs in intermediate mass halos of M 200 ~ 10^“^ hi©- 
Several groups have recently performed systematic 
searches through wide-area imaging surveys to locate 
such lenses with intermediate Eins tein radii — 2"5 — 
6 

Stark et al 

et al.| 20 d^ 


These include GASSOWARY (iBelokurov et al. |2009 

li 


et al. 2009 


2010 


the Sloan Bright Arcs Survey (|Dieh 
and the SL 2 S-JnTC^ 


sample (jLimousin et al. 2009| |More et al.| |2012|). We 
refer to these as group-scale lerises in r efere nce to their 


Einstein radii and halo masses (Section 3.3) that lie be¬ 
tween galaxy- and cluster-scale lenses. This term does 
not imply that galaxy-scale lenses residing in lower-mass 
halos are not found in the enriched environments t ypical 


of massiye galaxies, as many studies have shown (Kee- 


ton et al. 1120001 

IWilliams et al.| 

20061 lAuger et al. 1120071 

h'assnacht et al.||2008 

Treu et a 

ri[2ooir —— 

h'ollow-up studies oJ 

these group-scale lenses have lo- 


cuse d on their halo masse s and bulk mass-to- light ra¬ 
il^, scaling 

15m 


tio I 

Than 

avur et al. 

2010 

relations ( 

h'oex et al. 



the mass-concentrationTelation (iVLGR 


2013 Verdugo et al. ' 5D14), and 


2011 Wiesner et al.|[20T2} [Auger et al.|^13[ |Eoex et al. 


o et al. 


2014])! However, very few of these group-scale lenses 

have been studied using the lensing and stellar dynamics 
approach to measure the mass distribution within i?e ■ 
McKean et al. (2010) presented a detailed analysis of 
one radio-selected group-scale lens. Resolved stellar kine¬ 
matic data are needed to separate the stellar and dark 


been studie d in this way (Spiniello et al. 2011 Deason 

erar][20l^ . 


In this paper we present new observations of a sample 
of 8 group-scale lenses using the DEIMOS spectrograph 
at the Keck II telescope. The data allow us to the mea¬ 
sure the radial stellar velocity dispersion profile of the 
brightest group galaxy (BGG) and to estimate the halo 
mass based on the kinematics of the satellite galaxies. In 
conjunction with strong lensing, these data provide mass 
measures at several widely separated radii. We combine 
the new sam ple with earlier da,ta co llected by jSpinielloj 
et al. (2011) and Deason et al. (2013) to create a sample 
of 10 well-studied lenses. 'I'he resulting sample hlls in a 
long-standing gap in halo mass distribution of similarly 
analyzed lenses, enabling us to explore trends in the mass 
profiles of ETGs over halo masses of 10^^ Mq — 10^^ Mq 
using homogeneous techniques and data. 

An outline of the paper follows. The reader interested 
in only the results and not the methodology may wish to 
begin in Section 7. We introduce the sample in Section 2. 
Sections 3-5 describe the observations and the associated 
constraints from kinematics (§3), lensing (§4), and stellar 
population synthesis (§5). In Section 6 , we describe our 
procedure for inferring mass models from these various 
data sets. In Section 7, we begin presenting our results 
with the MGR for our group-scale lenses. In Section 8 , 
we move to smaller scales and examine the total density 
profile within Re, combining our new group sample with 
earlier work on galaxy- and cluster-scale lenses to study 
trends over a factor of ~ 60 in halo mass. In Section 9 
we examine these trends using a set of CDM-motivated 
halo occupation models. We investigate trends in the 
DM density profile within Re using both these halo oc¬ 
cupation models and direct mass modeling based on our 
resolved stellar kinematic data. Finally, in Section 10 we 
discuss the physical implications of our findings for mod¬ 
els of ETG formation and the effect of baryons on their 
structure of their DM halos, before summarizing our key 
results in Section 11. 

Throughout we use a ACDM cosmology with = 
0.3, = 0.7, and h = 0.7. All magnitudes are e xpressed 

in the AB system. Stellar masses are based on a jSalpeterj 
(19551 IMF over 0.1 — 100M©, and halo masses M 200 are 
dehned relative to the critical density. We adopt a cos- 
mological baryon fraction — 0.15 where necessary 

( Planck Gollaboration et al.||2014 ). 


2. GROUP LENS SAMPLE 

Our sample consists of 10 lenses selected to have Ein¬ 
stein radii 0Ein in the range « 275 — 5", whose basic char¬ 
acteristics and discovery references are listed in Table 
For CSWA163 and GSWAl we i ncorporated data pu^ 
lished byjDeason et al.|(|2013D andjSpiniello et al.|(|2011|), 
while our analyses ot the remaining 8 lenses are based on 
new observations. Figure compares the present sam¬ 
ple to the galaxy- and cluster-scale lenses that formed the 
basis of earlier joint analyses of strong lensing and stellar 
kinematics. The group-scale lenses have Einstein radii, 
stellar masses, and effec tive radii that bridge these earlier 
samples. In Section jTSj we show the average halo mass of 
our sample is M 200 = 10^"^ ° Mq, indicative of a group or 
small cluster betwe en the halo masses o f the galaxy-scale 
(M 200 — Mq; Gavazzi et al. 2007) and cluster-scale 
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aciq [km/sl 

200 300 500 700 900 



log ^Ein (Dds/D3=0.7) 

Fig. 1. — Characteristics of ETGs for which strong lensing and 
stellar dynamics have been combined to measure the density profile 
within ~ Re^ as a function of angular Einstein radius ^Ein- The 
present group-scale s ample bridges the large s amples of galaxy- 
scale len ses (SLAGS, |Auger et aL||2009| SL2S, |Sonnenfeld et al.| 
pOlSalb] ! with the massive clusters analyzed by INi^a, All 

yEin liave been rescaled to a common D^s/Ds = 0.7, and stellar 
masses M* have been homogenized to a Salpeter IMF. The top 
axis shows the velocity dispersion for a singular isothermal sphere 
with the indicated ^Ein- 

(M 200 = Mq; N13a) lenses in Figurej^ 

We required that our selected lenses be dominated by a 
BGG, since it is critical for our dynamical analysis that 
this galaxy be centrally located within the group-scale 
halo. Figure shows that the strong lensing region is 
dominated by a single deflector with the partial exception 
of CSWA141, which has a nearby satellite 1 mag fainter. 
On larger scales, we find that the central group galaxies 
are 0.9-1.7 mag brighter than the second-rank candidate 
group members]^ 

3. SPECTROSCOPIC DATA 

Here we present spectroscopic observations of the 
BGGs and candidate group members. Our goals are (1) 
to measure the internal kinematics of the BGG through a 
radial velocity dispersion profile, and (2) to identify other 
group members whose velocities will constrain the halo 
mass. We also search for other line-of-sight s truc tures 
that will inform our lensing analysis in Section |4.2[ 

3.1. Observations and Reduction 


^ We define these as the second-brightest galaxy within 500 kpc 
that has a pho tometric redshift in the Sloan Digital Sky Survey 
DRIO catalog ||Ahn et al.||20i^ within 0.1 of the lens galaxy. 
J09413 lies outside t he SP SS footprint, so we instead use the images 
described in Section|3.1|and compare to red sequence members. 


We designed slit masks for the DEIMOS spectrograph 
(Faber et al. 2003) at_Ure Keck II telescope for the 8 
groups listed in I'able targeting both the BGG and 
candidate satellites. Two masks were designed for 4 
of the groups, with the satellite targets switched while 
the BGG slit was fixed on both masks. Three other 
groups were observed with a single slit mask, while the 
CSWA165 BGG was observed with a long slit. The 
600 mm“^ grating was used in combination with the 
GG455 or GG495 blocking filters and a 1" slit. Gan- 
didate satellites were drawn from Sloan Digital Sky Sur¬ 
vey (SDSS) catalogs (Ahn et al. 2014), prioritizing red 
sequence members, followed by galaxies with photomet¬ 
ric redshifts consistent with being members, and finally 
other bright galaxies in the field. Observations of the 
8 groups were conducted over 4 nights in 2013 in clear 
conditions and seeing of O'.'? — 0'.'9. Total exposure times 
on the BGGs were 3-4 hr. 

The data were reduced using the spe c2d pipeline 
(Cooper et al. 2012 Newman et al. 2013). Since the 
default sky subtraction routines are not appropriate for 
the extended BGGs, we adapted them to accommo¬ 
date a more generous mask of the galaxy light. Red- 
shifts of the satellite candidates were measured by cross- 
correlating wit h absorption- and emis sion-line templates 
using rvsao (Kurtz & Mink 1998) for 647 galaxies, 
whose velocities relative to the BGG were tabulated as 
Au = c{z - 2:bgg)/(1 + -2:bgg)- 


3.2. Lens Environments 

We first use our redshift surveys to probe the large- 
scale environment of the groups and determine whether 
there are any line of sight stru cture s relevant for our 
strong lensing analysis (Section |4.2[ ), recognizing that 
our survey may be incomplete as it was intentionally bi¬ 
ased toward group members. Nevertheless, structures 
are found in several cases. Figure shows the redshift 
distribution around each of the 7 lenses that were ob¬ 
served in multi-slit modej^ CSWA165, CSWA107, and 
CSWA141 show no sign of additional structures in the 
field. J09413 is likewise dominated by the lensing group, 
with only a mild secondary peak that is located far from 
the lens (filled histograms show galaxies within 1 Mpc of 
the BGG). CSWA6 and CSWA7 each overlay comparably 
rich structures. EOCL is the most complex system, with 
three r edsh ift peaks within 4000 km s“^ of the lens. In 
Section [4^ we use these results to judge the contribution 
of external structures to our lens models. 


3.3. Satellite Kinematics and Halo Masses 

The kinematics of the satellite galaxies of the group 
provide a measure of the mass on scales extending to 
virial radius. The first step in such an analysis is to 
identify the group members. We adopted an iterative 
cut in phase space that rejects galaxies with |u —wbggI > 
3cr(i?), where each velocity is compared to the local ve¬ 
locity dispersion a(R) appropriate to its projected group¬ 
centric radius R. For this purpose we set the shape of the 
<j{R) profile using a fiducial NFW halo with concentra¬ 
tion c = 4 and isotropic orbits, while the overall velocity 

® Since the limiting magnitude of the SDSS catalogs reach fainter 
galaxies in lower redshift systems, this figure should not be used 
to gauge the relative richness of the groups. 
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Fig. 2.— Le/t panels: Images of the inner regions of each lens from the sources described in Section [4.1| Rulers in the bottom-left corner 
have a len gth o f 3”; note the scale varies among images. Right panels: Images generated from the best-titting analytic lens models described 
in Section |4.2| The solid line shows the outer critical curve. For the case of CSWA7, the positions of the compact multiple images (green 
circles) are used to constrain the lens model rather than the pixel-level data; this accounts for the lack of an arc image in its right panel. 
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TABLE 1 
Lens Sample 


Name 

R.A. 

Dec. 

zl 

zs 

®Ein 

References 

CSWA107 

11:47:23.30 

+33:31:53.6 

0.212 

1.205 

2'.'52 

S13 

CSWA141 

08:46:47.46 

+04:46:05.1 

0.241 

1.425 

3'.'15 

S13 

CSWA164 

02:32:49.87 

-03:23:26.0 

0.450 

2.518 

3'.'68 

S13 

CSWA165 

01:05:19.65 

+01:44:56.4 

0.361 

2.127 

4'.'33 

S13 

CSWA6 (The Clone) 

12:06:02.09 

+51:42:29.5 

0.433 

2.00 

4'.'36 

L09, S13 

CSWA7 

11:37:40.06 

+49:36:35.5 

0.448 

1.411 

2'.'73 

K09, S13 

8 O’clock Arc (EOCL) 

00:22:40.91 

+ 14:31:10.4 

0.380 

2.73 

3'.'29 

A07 

J09413-1100 

09:41:34.7 

-11:00:54.3 

0.385 


4'.'04 

Li09 

CSWA163 

21:58:43.67 

+02:57:30.2 

0.287 

2.081 

3'.'49 

S13, D13 

Cosmic Horseshoe (CSWAl) 

11:48:33.14 

+ 19:30:03.1 

0.444 

2.379 

5'.'08 

B07, D08, Sll 


Note. — ^Ein is the Einstein radius as mea sured in Sect i on 14.2 [ zi, a nd z<^ are t he len s and 
source redshifts , respe ctively. References: A 07: Al lam et al.|(^U7|l, D OS: i Dye et al.||l2008 , L09 : 
iLin e t a'n<|2009L Li09: iLim ousi n et al.|l|2009ll, K Jij: llAubo et ahld^UOQL Sll: |bpinielio et ai.| 2011|l, 
JJlIj: [jJ^son et al.| ( |20i;j| |, aid: [iStark et alT[ l |2013| l 
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Fig. 3. — Redshift distribution in the vicinity of each lens derived from our spectroscopic survey. All galaxies in the DEIMOS field of 
view are plotted. Dark colored histograms show the subset within R < 1 Mpc of the BGG^while the light blue histograms show the subset 
of these that are also identified as group members by the R — v criterion shown in Figure!^ 
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Fig. 4.— Distribution of galaxy velocities relative to each BGG, v = c{z — 2bgg)/(1 + 2:bgg)i a-s a function of projected separa tion. 
Dashed curves show the region within which group members (filled circles) are selected using the iterative procedure described in Section [3.3| 
These members are used to calculate the satellite velocity dispersion cr, which in turn is used to estimate the halo mass M2oo- 
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TABLE 2 

DEIMOS Spectroscopic Observations Log 


Name 

Dates 

£exp 

(min.) 

Masks 

Redshifts 

P.A. 

(deg) 

CSWA107 

2013 Nov 27-28 

180 

2 

102 

83 

CSWA141 

2013 Nov 27-28 

210 

2 

93 

37 

CSWA164 

2013 Feb 10-11, 
2013 Nov 28 

242 

Oa 


-74 

CSWA165 

2013 Nov 27-28 

196 

2 

78 

-71 

CSWA6 

2013 Feb 11 

180 

1 

60 

25 

CSWA7 

2013 Feb 10 

230 

2 

102 

75 

EOCL 

2013 Nov 27 

210 

2 

91 

12 

J09413 

2013 Feb 10-11 

240 

1 

62 

-30 


et al. (2013) adapted the tracer mass estimator (TME) 
formulated by Evans et aL (|2003 1, which they calibrated 
to iV-body siinuiations. |Zhang et ai.| (|2011|) measured 
an empirical scaling relation between a and r^oo based 
on X-ray measurements}^ Using the TME or Zhang et 
al. mass estimators shifts the halo masses systematically 
by —0.05 dex and -1-0.08 dex, respectively. We con¬ 
sider 0.08 dex as a reasonable estimate of the system¬ 
atic uncertainty in the mass scale, which we add to the 
random errors in a. The final halo masses and their 
uncertainties are listed in Table and span the range 
logMaoo/Mo = 13.7-14.5. 


Note. — Exposure time £exp refers to the total integration on 
the lens galaxy. 

^ CSWA164 was observed using a long slit. 


TABLE 3 

Satellite Kinematics 


Name 

-^gal 

cr 

(km s“^) 

(v) 

(km s“^) 

log M200I Mq 

CSWA107 

45 

552 ± 59 

38 ±91 

14.18 ±0.18 

CSWA141 

48 

374 ± 39 

45 ±64 

13.76 ±0.17 

CSWA165 

26 

362 ± 52 

82 ±61 

13.78 ± 0.21 

CSWA6 

12 

456 ± 101 

124 ± 175 

14.07 ± 0.28 

CSWA7 

15 

384 ± 75 

-115± 111 

13.89 ±0.25 

EOCL 

14 

319 ± 65 

3± 124 

13.65 ±0.26 

J09413 

22 

655 ± 104 

144 ± 118 

14.45 ±0.22 

CSWA163 

22 

654 ± 103 

-35 ± 162 

14.40 ± 0.22 


Note. — Agai is the number of spectroscopically identified 
group members, (v) is the mean velocity of the satellites with 
respect to the BGG. 


scale is then set by matching the aperture velocity dis- 
persion within the v irial radius to the measured value. 


Mamon et al. (2013) advocate a similar cleaning proce- 
dure and provide useful analytic approxi mations (see also 


Katgert ^ al.|[2004 Biviano et ai.||2006 ). 
h'igurel^shows the R—v plane tor the eight groups with 


mea sure Jsatelli te kinematics (seve n introduced in Sec¬ 
tion ITT plus the Deason et al.||2013 data for CSWA163). 


The curves show the '6a threshold for selecting group 
members identified with filled symbols. The effective¬ 
ness of this procedure is demonstrated by the rejection of 
some galaxies that might have been included by a simple 
velocity cut (in particular, several in CSWA7 and EOCL) 
which are found to be spatially coherent substructures. 
We then calculate the line-of-sight velocity dispersion a 
of the members using a simple standard deviation with 
an uncertainty estimated from Monte Carlo simulations. 
Table 1^ shows that a ranges from 319 to 655 km s“^, 
with a median of 455 km s“^. Furthermore, the BGGs 
are consistent with being at rest with respect to their 
satellites as expected if they are centrally located in their 
halos. 

To estimate the halo mass, we use the sca ling re¬ 
lation betwe en a and M 2 Qn determined by |Munari 


et al. (2013) in simulations^ ^og h(z)M 2 oo = 13.98 -|- 
2.75 log cr/(500 km s“^). We compared these mas ses with 
those obtained from two alternate approaches. |Deason| 


3.4. Stellar Kinematics of the Central Galaxy 

Spatially resolved spectra of the BGGs were extracted 
and analyzed to derive a radial profile of the projected 
stellar velocity dispersion a\os{R) following the proce¬ 
dures described in N13a. Briefly, extraction bins were 
constructed to ensure a minimum signal-to-noise ratio 
of 15 A-i in the rest frame in the wavelength range 
4150 — 4950 A around the G band. Kinematics were 


measured using ppxf (Cappellari & Emsellem 2004). Op¬ 
timal stellar templates were constructed from a linear 
combination of spectra of G and K giants with metallic- 


ities near solar drawn from the MILES library (Sanchez- 
Blazquez et al. 2006). The templates were redshitted 


convolved with a Gaussian (taking into account the in¬ 
strumental resolution tr = 78 km s“^), and added to 
a polynomial to filter the continuum, following well- 
established procedures. As a typical example. Figure 
demonstrates the high quality of the resulting fits for 
GSWA7. We estimate a 5% systematic uncertainty in a 
by varying the fitted wavelength region, the polynomial 
order, and the template library. Since the uncertainty 
is highly correlated amongst the spatial bins, we do not 
add this in quadrature to the random errors, but include 
a calibration factor with a Gaussian prior in our mass 
models (Section}^. Figure shows the derived veloc¬ 
ity and velocity dispersion profiles. Rotational support 
is negligible or absent in every case, so we ignore it in 
our dynamical modeling. Our stellar velocity dispersion 
measurements are listed in Table ID 


4. LENS MODELS AND GALAXY SURFACE PHOTOMETRY 

We now turn to our method for analyzing the strong 
lensing. After introducing the imaging data, we describe 
our technique for fitting the images at the pixel level us¬ 
ing analytic models for the mass and light distributions 
of the lens and source. This allows a precise measure¬ 
ment of the Einstein radius. At the same time, we obtain 
multicolor surface brightness profiles of the BGGs, a key 
ingredient for our mass modeling in Section 

4.1. Imaging Data 

GSWA6, GSWA7, and EOGL have been imaged by 
WFPG2 onboard the Hubble Space Telescope (HST) 
through the F450W, F606W, and F814W filters (pro¬ 
gram IDs 11167 and 11974, P.I. Allam). J09413 was 
imaged using HST/ACS through the F475W, F606W, 


_^_Wetake^h«caliLration_intheiT_^Llela££ro£riate_^o^al^w 

tra^erainsimuJatinn^withAG^teedbacUVa^inKth^ra^erand 

pliysic'sni3TgesTEe''caIi^ra^ori^^'^^LT58aSc^^^”””” 


® Here we use their BCES bisector fit to the whole sample and 
convert M200 = l.SSMsoo; this conversion is exact for NEW halos 
with c = 5 and depends weakly on concentration. 
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Rest wavelength [A] 


4200 4300 4400 4500 4600 4700 4800 4900 



Observed wavelength [A] 


Fig. 5. — Example of our resolved spectroscopy of BGGs in the case of GSWA7, showing the data (gray) and the convolved stellar 
template (red) used to extract the velocity dispersions in various spatial bins indicated in the figure. These bins are symmetric about 
the center of the galaxy. Spectra are smoothed with a 5 A boxcar, normalized to a median flux of unity, and offset vertically for clarity. 
Spectral regions with uncertain calibrations were excluded from the fit and are not drawn in the model fits. Errors in a are statistical and 
do not include the systematic uncertainty of 5% described in the text. 



Radius [kpc] 


Radius [kpc] 


Fig. 6. — Stellar kinematics of the BGGs extracted in spatial bins along the spectroscopic slit. Black crosses with error bars show the 
measured velocity dispersions (Jios, while the red solid and dashed lines indicate the velocity profiles and their ±lc7 uncertainties. Velocity 
profiles are shifted vertically so that the systemic velocity is at +350 km s“^. The gray band shows the crios profiles of the fitted mass 
models, introduced in Section]^ and enclose 68% of the posterior samples. 
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TABLE 4 

Velocity Dispersion Profiles of BGGs 


Name 

Bin limits I 

arcsec) 

cr (km s 

Name 

Bin limits ( 

arcsec) 

(7 (km s 

CSWA107 

-3.50 

-2.67 

269 ± 21 

CSWA6 

-1.48 

-0.89 

331 ± 21 


-2.67 

-2.07 

300 ± 20 


-0.89 

-0.30 

309 ± 10 


-2.07 

-1.48 

265 ± 13 


-0.30 

0.30 

299 ±8 


-1.48 

-0.89 

223 ±9 


0.30 

0.89 

308 ±9 


-0.89 

-0.30 

257 ± 6 


0.89 

1.48 

349 ± 19 


-0.30 

0.30 

254 ± 6 


1.48 

2.43 

406 ± 42 


0.30 

0.89 

247 ± 6 

CSWA7 

-2.07 

-1.48 

258 ± 23 


0.89 

1.48 

261 ± 10 


-1.48 

-0.89 

260 ± 15 


1.48 

2.07 

248 ± 15 


-0.89 

-0.30 

264 ± 7 


2.07 

2.67 

221 ± 17 


-0.30 

0.30 

258 ±6 

CSWA141 

-3.02 

-2.07 

414 ± 36 


0.30 

0.89 

259 ±8 


-2.07 

-1.48 

337 ± 27 


0.89 

1.48 

239 ± 20 


-1.48 

-0.89 

258 ±9 


1.48 

2.07 

269 ± 28 


-0.89 

-0.30 

247 ± 8 

EOCL 

-3.02 

-2.07 

302 ± 39 


-0.30 

0.30 

236 ± 5 


-2.07 

-1.48 

292 ± 33 


0.30 

0.89 

252 ± 7 


-1.48 

-0.89 

270 ± 17 


0.89 

1.48 

262 ± 11 


-0.89 

-0.30 

298 ± 7 


1.48 

2.07 

281 ± 26 


-0.30 

0.30 

309 ±8 


2.07 

3.02 

363 ± 34 


0.30 

0.89 

279 ±9 

CSWA164 

-1.60 

-0.89 

288 ± 37 


0.89 

1.48 

285 ± 26 


-0.89 

-0.30 

326 ± 17 

J09413 

-2.31 

-1.48 

291 ± 27 


-0.30 

0.30 

314 ± 13 


-1.48 

-0.89 

344 ± 17 


0.30 

0.89 

303 ± 17 


-0.89 

-0.30 

334 ± 10 


0.89 

1.60 

337 ± 45 


-0.30 

0.30 

351 ± 8 

CSWA165 

-2.19 

-1.48 

343 ± 32 


0.30 

0.89 

333 ±9 


-1.48 

-0.89 

287 ± 17 


0.89 

1.48 

326 ± 16 


-0.30 

0.30 

281 ± 9 


1.48 

2.31 

378 ± 38 


0.30 

0.89 

308 ± 9 






0.89 

1.48 

326 ± 17 






Note. — Velocity dispersions are measured in rectangular apertures defined by the 1" slit width and 
the bin limits along the slit, which are tabulated relative to the galax y center. Errors in a are s tatistical 
only and do not include the estimated 5% systematic uncertainty. See|Deason et al.|(|2013|l and|Spiniello| 
|et al.| | [2011[ |, respectively, for the CSWA163 and CSWAl data. 
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and F814W filters (ID 10876, P.I. Kneib). CSWAl was 
observed with HST /WFC3-UVIS through these same fil¬ 
ters (ID 11602, P.I. Allam). For CSWA163 we rely on gri 
images from the SDSS. As part of the Keck/DEIMOS ob¬ 
servations described in Section [3A| we imaged CSWA107, 
CSWA141, CSWA164, and CSWA165 through the B and 
R filters. Exposure times ranged from 4 to 12 minutes 
with seeing of (f!7 — I'.'O. Astrometric and photomet¬ 
ric solutions for the DEIMOS images were derived from 
stars in the SDSS catalog. 


4.2. Modeling the Lens Systems 

We fit analytic models of the mass and light distribu¬ 
tion to the multi-band data introduced above. By di¬ 
rectly fitting the image pixels in several filters simulta¬ 
neously, we naturally de-blend the lens and source galaxy 
light. Eigurel^shows the regions around each lens used to 
fit the modeb^which roughly encompass the radii where 
uncertainties in the background level have a minimal ef¬ 
fect. Light from the lens and source galaxies is modeled 
with a seven-parameter elliptical Sersic profile. For the 
BGG, we let the magnitude and Re vary among filters 
to allow for color gradients, but we fit a common Sersic 
index n, center {xo,yo), position angle (P.A.), and axis 
ratio q = b/a. For the background source, only the mag¬ 
nitude is allowed to vary among filters. 

The deflecting mass is modeled as a power law pro¬ 
file characterized by a slope 7 sl, where p oc and 

an Einstein radius 0Ein- Ellipticity is intr oduced into the 
surfa ce density following Schramm (1990 see also Keeton 
20011. The P.A., axis ratio, and center of the mass dis¬ 
tribution are generally not tied to those of the light. The 
exceptions are the naked-cusp configurations (GSWA107, 
CSWA7, J09413) and GSWA141 (see Appendix A) for 
which we found that the center and ellipticity cannot 
both be constrained. In these cases, we fixed the center 
of mass to that of the BGG. External shear with ampli¬ 
tude F and orientation Oy is also incorporated in the lens 
models where the data quality and image configuration 
can provide sufficiently useful constraints; this is the case 
except for CSWA107, GSWA141, and GSWA163. 

In addition to the main deflector, we model the de¬ 
flection from satellite galaxies in CSWA165, GSWA107, 
CSWA141, and GSWA6 where the perturbing galaxies 
are clearly visible in Figure!^ The satellite galaxy light 
was again modeled with a sersic profile, while the mass 
was treated as a singular isothermal ellipsoid (SIE) whose 
center, axis ratio, and P.A. were fixed to those of the 
light. This leaves a single free parameter, a. We place 
a Gaussian prior on a us ing the [Faber fc Jacks on] (fl976 1 
relation, as measured by Bernardi et al. 1 2UU3). similar 
procedures are commonly used m cluster fens modeling 
(see N13a and references therein). 

We define the Einstein radius 0Ein such that 
Atot(^*Ein) = 1, where Ktot is the azimuthally averaged 
mean convergence profile. We then define Kgroup as the 
convergence within 0Ein associated with the main deflec¬ 
tor, i.e., the BGG and group-scale halo. This differs 
slightly from unity for the four lenses with perturbing 


satellites included in the model, k. 


:rou p 


then serves as 


input to the mass models in Section |6P| We describe our 


® Although the mass models in Section subdivide the main 


treatment of external convergence /text along with our 
mass modeling procedure in Section [6.1| 

To fit the lens m odels, the paramet er space is explored 
using MultiNest (Feroz et al. 20091, a Markov Ghain 
Monte Garlo engine. For a given set of parameters, we 
generate galaxy images in the lens plane and ray-trace 
the source galaxy through the mass distribution. Images 
are generated for each observed filter and are convolved 
by the relevant point spread function (PSF). These model 
images are compared to the data to compute a likelihood 
L oc exp(—|x^). Figure shows that the best-fitting 
models provide acceptable fits. Although simple ana¬ 
lytic models cannot be expected to trace the detailed 
source structure in all cases (e.g., GSWAl), they are ad¬ 
equate for measuring 0Ein- Because its multiple images 
are virtually unresolved, CSWA7 was analyzed using a 
different technique that incorporates only the astromet¬ 
ric positions of these images (see Appendix A). 

With high-quality imaging data, the formal statisti¬ 
cal errors on the mass model are very small; system¬ 
atic differences are much more important. To assess 
these, we first modeled th e lenses using a non-si ngular 
isothermal ellipsoid (NIE; Kormann et al. 1994). We 
also tested models excluding external shear! finally, we 
fit bo th NIE and power law models using the glafic 
code ( Oguri]20l0 ) and the positions of conjugate images 
as constraints, rather than the pixel-level data. By com¬ 
paring these different methods, we find that 7%hiup varies 
by less than about 0.05, which we take as a fiducial un¬ 
certainty. Higher uncertainties adopted in a few cases are 
discussed in Appendix A, where circumstances particu¬ 
larly to individual systems are reviewed and comparisons 
are made to models in previous publications. 

The resulting lens model parameters are listed in Ta¬ 
ble We generally find good alignment between mass 
and fight: (1) spatial offsets are < O'.'! for all systems ex¬ 
cept GSWA6, (2) the axis ratios q agree on average, with 
a scatter of 0.1, (3) the position angles agree with a scat¬ 
ter of 16°. These comparisons support our assumption 
that the BGGs are centrally located within the group- 
scale halos. 


4.3. BGG Surface Photometry 


As described in Section [4~1] we fit Sersic profiles to the 
BGG surface photometry in several filters as part of the 
lens modeling procedureH These profiles were then in¬ 
terpolated to rest-frame .o and V filters. The B filter en¬ 
compasses the G band region where the BGG kinematics 
were measured, while the V filter is the reddest gener¬ 
ally available and should better trace the stellar mass; 
both are needed for our mass modeling. For each lens, 
we fit a linear relation to the fc-correction derived from 


Bruzual & Chariot (2003 BC03) simple stellar popula- 
tion models as a tunction of the observed color, using 
the observed filter pair nearest to the redshifted B or 
V bands. After removing Galactic extinction following 


deflector into its stars and DM halo, the ^Ein derived from lensing 
using single power-law models is still valid, sinc e is known to 
be nearly independent of the mass profile (e.g., |Rusin et al.|2003| 
and tests in Section 4.2). 

^ For CSWA164, we found it necessary to perform a separate 
fit to the BGG light alone after masking the Einstein ring, as it is 
significantly blended with the BGG in the ground-based imaging. 
The Sersic indices were also allowed to vary between Alters. 
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TABLE 5 

Strong Lensing Constraints 


Name 

^Ein 

^group (f^Ein ) 

7SL 

Qinass 

i^A-mass 

(ARA, ADec) 

Text 

Or 

GSWA107 

2'.'52 

0.96 ±0.10 

1.18 

0.73 

-89.6 




GSWA141 

3'.T5 

0.91 ± 0.05 

1.23 

(fixed to BGG) 




GSWA164 

3'.'68 

1 ± 0.05 

1.51 

0.88 

-18.5 

(0i'01,0i'00) 

0.022 

80.8 

GSWA165 

4'.'33 

0.97 ± 0.05 

1.73 

0.77 

-72.3 

(0i'02,0i'01) 

0.056 

-66.6 

GSWA6 

4'.'36 

0.79 ±0.10 

1.79 

0.86 

-60.2 

(-0727,0730) 

0.012 

-46.4 

GSWA7 

2'.'73 

1 ± 0.05 

1.85 

0.67 

61.6 


0.149 

71.4 

EOCL 

3'.'29 

1 ± 0.05 

1.91 

0.74 

14.9 

(0i'01,0i'07) 

0.072 

10.8 

J09413 

4'.'04 

1 ± 0.07 

1.59 

0.54 

-14.2 


0.087 

57.0 

GSWA163 

3'.'49 

1 ± 0.05 

1.74 

0.71 

87.9 

(0i'02,0i'02) 



GSWAl 

5'.'08 

1 ± 0.05 

1.66 

0.90 

-52.9 

(Oi'03,0706) 

0.022 

-5.7 


Note. — Kgroup is the azimuthally averaged mean convergence of the main deflector 
measured within the Einstein radius fiEin- (For a single deflector, this is unity by definition, 
but differences arise when perturbing galaxies contribute convergence.) Offsets (ARA, ADec) 
give the center of mass relative to that of light; where omitted, the center of mass is fixed. 
In models with external shear. Text and specify its amplitude and orientation (east of 
north). 


TABLE 6 

Sersic Profile Fits to BGGs and Mass-to-Light Ratios from SPS Modeling 


Name 

q 

RA. 

(deg) 

n 

Re,B 

(kpc) 

0 

o 

Re,V 

(kpc) 

0 

o 

tSPS 
^ V,0 

A^©/-^0 

(T|fS) 

CSWA107 

0.60 

87.6 

4.01 

31.6 

2.46 

25.8 

2.40 

5.1 

3.7 

CSWA141 

0.79 

54.3 

5.84 

56.0 

3.31 

40.3 

2.90 

5.1 

3.7 

CSWA164 

0.89 

-26.9 

3.15 

13.4 

2.50 

10.8 

2.39 

3.2 

3.0 

CSWA165 

0.78 

-74.6 

6.10 

42.1 

3.32 

29.6 

3.29 

4.9 

3.8 

CSWA6 

0.85 

-86.1 

6.51 

36.8 

3.67 

30.0 

3.66 

4.6 

3.7 

CSWA7 

0.65 

63.8 

6.07 

29.9 

2.49 

24.2 

2.47 

4.5 

3.7 

EOGL 

0.74 

10.3 

5.76 

31.3 

2.74 

25.7 

2.79 

4.7 

3.8 

J09413 

0.68 

-6.5 

6.89 

53.4 

4.00 

46.1 

4.60 

5.2 

4.0 

CSWA163‘‘ 

0.81 

-84.3 

4t 

21.3 

2.00 

17.6 

2.04 

4.7 

3.8 

CSWAl 

0.90 

-15.3 

5.79 

26.7 

3.14 

20.9 

3.35 

3.9 

3.3 


Note. — B and V refer to the rest frame. Radii are circularized. T^^q® specifies Mi,ILy 
measured withi n the central spectroscopic aperture (If'O X 0”6) based on SPS models and a 
|Salpeter| | |1955^ IMF. The li^t- weighted mean Mf/Ly within the E-band effective radius is 
(T^*^^). See Sections 4.3 and nil for a discussion of the uncertainties. 

^ The Sersic index wasnxed ter CSWA163 due to the poorer quality of the SDSS imaging. 


Schlafly & Finkbeiner (2011|), we apply this radially de- 
pendent fc-correction and shift to the rest frame to obtain 
a surface brightness profile in units of Lq kpc“^. This 
is then fit with a Sersic model having a free Re and to¬ 
tal luminosity, but with n fixed to the value measured 
in the observed-frame fits. Table [6] lists the rest-frame 
surface brightness profiles for each BGG. Errors in the 
Sersic parameters are highly correlated. The most rele¬ 
vant measure of uncertainty for our analysis is the am¬ 
plitude of systematic deviations between the observed 
and model surface brightness profiles, which is typically 
< 0.15 mag arcsec“^ within i? < 8". 


5. STELLAR POPULATION SYNTHESIS AND RADIAL 
GRADIENTS 


One goal of our analysis is to compare the stellar mass 
obtained from lensing and dynamics with that estimated 
with stellar population synthesis (SPS) models. In this 
section we analyze spectroscopic and photometric obser¬ 
vations of the BGGs using SPS models to constrain the 
stellar mass-to-light ratio = JVR/Ly. 


As a first step, we estimate a central value The 

DEIMOS spectrum extracted from the central 1"0 x 0"6 


of each BGG is fit using the pyspecfit code (Newman 


et al 


1[M4l ) 


and a suite of BC03 models based on a fidu- 
cial Salpeter IMF. The models follow exponentially de¬ 


clining star formation histories with uniform pri¬ 

ors on 0 < logage/Gyr < 1, 7 < logr/yr < 10, and 
0.01 < Z < 0.04, where Z is the metallicity. The red- 
shift and velocity dispersion were also fitted simultane¬ 
ously, and a 12th order multiplicative polynomial was 
used to filter the continuum. Figure shows the result¬ 
ing fits and the uniformity of the I^G spectra. This 
leads to a narrow range of Ty(f estimates listed in Ta¬ 
ble The uncertainties are dominated by systematics 
discussed below. 

A common approximation in lensing and dynamical 
studies of ellipticals is that stellar mass follows the optical 
light profile, with a radially invariant M* jL. However, 
Figure [^demonstrates that all of the BGGs in our sample 
show negative color gradients, which implies that A'R/L 
declines with increasing radius. To quantify this decline, 
we construct (7, B, and V band rest-frame Sersic profile 
fits as described in Section |T3| For each BGG, we choose 
the rest-frame color most closely matching one of the 
observed pairs, which ensures that only a small color 
interpolation is needed for most systems (CSWA164 and 
CSWA165 are exceptions). The color gradient is then 
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Fig. 7.— Spectroscopy of the centers of the 8 BGGs observed 
with DEIMOS, extracted in a l70 x 076 aperture (gray) and fitted 
with SPS models (red) to estimate as described in Section!^ 

The spectra are highly uniform. Regions contaminated by bright 
sky lines or uncertain relative flux calibration are not displayed. 


fitted by a linear slope Vcolor = dcolor/dlogi? over the 
interval 0"3 < R < 5", which encompasses our leasing 
and dynamics constraints. 

Metallicity gradients are generally found to be more 


Tamura et al.||2000l Mehlert et al. 

2003 Wu et al. 

Kawle et al.|2fllUt|'I'ortora et al.|2U] 

(lllGreene et al. 


color to M^/Ly gradients, we use relations derived from 
BC03 simple stellar population models with a fixed age of 
7 Gyr and varying metallicity: dlogTv/d{B — V) = 1.2, 
d\ogTy/d{U—B) — 0.72, and dlogTy/d{U—V) = 0.45. 
(Note that with more extensive photometry extending 
into the the near-infrared, we could address this de gen¬ 
eracy direct ly in these objects.) For comparison, the 


Bell 


et al. (2003) relation — which is based on trends in age, 
dust, and rnetallicity found empirically from multi-band 
spectral energy distribution fitting of the overall galaxy 
population — gives a very similar cilogT\//(i(i3 — V) = 

1.305. 


The resulting gradients VTy = d\og{M^/Ly)/dlogR 
are listed in FigureThe median is VTy = —0.15 with 
a dispersion of only 0.03j^ In our fiducial mass models, 
we therefore place a Gaussian prior on VTy with this 


® Here we have excluded CSWA164 and CSWA165, since their 
rest-frame colors required a significant extrapolation, as well as 
CSWA163 due to its lower data quality, but these exclusions turn 
out to have a minimal effect on the median. 


mean and dispersion. The radial variation is then 

, R \ VTv 

xSPS^n\ _ -wSPS ( ^ 

1 y \R) — 1 y,o 1 


( 1 ) 


since the mean radius in the central DEIMOS aperture 
is i? = O'.'S. Using Equation we can also compute the 
light-weighted mean wit hin the U-band effective 

radius, which we use in Section [973] to test the effect on 
our results of neglecting M */L gradientsj^ 

There are clearly significant systematic uncertainties 
in both the zeropoint of the estimates and their 

trend with radius. Concerning the former, we compared 
our measurements to t hose derived from the SP SS griz 
colors using kcorrect (|Blanton & Roweis||2007|), shifted 
to a Salpeter IMF. Since the colors are measured in large 
apertures, we compare the results to our (Ty^®). This 
approach uses different data and models and so provides 
a useful route to estimate uncertainties. We find an rms 
difference of 0.05 dex between the methods, which we 
take as a fiducial random uncertainty; the systematic 
shift between the methods is only 0.03 dex. 

Concerning the M^/L gradient, we have estimated a 
lower limit by supposing that the color gradient is domi¬ 
nated by metallicity variations: gradients in age or dust 
would give larger M^/L variations for a given color dif¬ 
ference. However, the color gradients in our sample 
are larger than the typical values se en in other stud- 
ies of elliptical galaxies. For example, Wu et al. (20051 
find a mean V(i? — U) = —0.05 ± O.OI, compared to 
a mean V{B — V) = —0.13 for the 6 g alaxies with 
that color p lotted in Figure Similarly, [Kuntschner 
et al.| (|2010|) measure a typicS spectroscopic metallic¬ 
ity gradient of d\ogZ/d\ogR = —0.25 corresponding to 
V(H -V) = -0.07 and VTy = -0.08 in the BC03 
models, which are about half our inferred values. It is 
not obvious why our sample would show much stronger 
stellar population gradients, although it may relate to 
differences in the assembly histories (e.g., properties of 
the cannibalized satellites) of central galaxies in 10^"^ Mq 
halos relative to typical “field” galaxies. 

In our mass models, we use Tyg® to inform broad priors 
on the stellar mass-to-light ratio, and we vary the radial 
gradient VTy within the range desc ribed above to test 
the sensitivity of our results (Section 9.3). 


6 . MODELING THE MASS DISTRIBUTION OF 
GROUP-SCALE LENSES 

With the observational constraints in hand, we now 
outline our method for inferring the mass distributions of 
the group-scale lenses. The key ingredients—the Einstein 
radius 0Ein, the velocity dispersion profile <j[R) of the 
BGG, and the halo mass M 200 of each lens—are used 
to constrain two-component mass models, consisting of 
a group-scale DM halo and the stellar mass of the BGG. 

The BGG is modeled using the Sersic fits presented in 
Section [T^ The stellar mass profile follows the U-band 
light multiplied by the radially dependent given by 


® For CSWAl and CSWA163, since we do not have access to the 
resolved spectra to measure we instead used pyspecfit to 

fit the SDSS griz photometry measured in a 20 kpc aperture. In 
conjunction with the mean gradient (VTy), we can then solve for 
the central Ty^Q® that reproduces this aperture measurement. 
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Fig. 8. — Color gradients in our sample of BGGs. Each panel shows the radial color variation as measured in directly in the images (blue 
and red diamonds) and in the Sersic model fits (dotted lines). Other galaxies, including the lensed object, are masked. Dashed lines show 
a rest-frame color obtained using the radially varying fc-corrections described in the text. The slopes of the color (Vcolor = dcolor/dlogit) 
and Ty (VTy = d log Ty/d log it) gradients are given in lower-right corner of each panel. 


Equation The luminosity profile at the wavelengths 
where kinematics are measured is also required to solve 
the Jeans equations, and for this we use the rest-frame 
B-band profile described in Section [473| The DM halo is 
modeled as a generalized NEW (gNFW) profile 


PBuir) = 


Po 


(r/rs)/5(l-|-r/rs)3-/3’ 


( 2 ) 


where the asymptotic inner slope /3 = 1 in the case of an 
NEW model. 

We compute the kinematic observabl es using the Jean s 
anisotropic modeling (JAM) routines ( Cappell^|2008 ), 
which operate quickly on oblate, axisymmetric mass dis- 
tributions. This allows us to move beyond spherically 
symmetric models often used in earlier work, including 
N13, and to test the effects of that assumption on our 
results. JAM requires that the surface density of mass 
and tracers be expressed as a multi-Gaussian expansion 
(MGE). Eor the Sersic and gNEW profiles, we therefore 
determined the MGE coefficients as polynomial functions 
of n and (3, thus allowing the MGE to be quickly con¬ 
structed for a given model 

The observed projected axis ratio q and the incli¬ 
nation angle i determine the intrinsic axis ratio q' = 
— cos^i /sini of the three-dimensional mass distri¬ 
bution. We take a uniform prior 0 < cosz < coszmin, 


For the gNFW profile, 18 Gaussian components were fit to 
pir/'f's) for many values of /3. The amplitudes of these components 
were then fit by 12th order polynomials in /3 while the widths 
Gi of the components were kept fixed. Over the range 10“^ < 
r/vs < 10 and 0 < /5 < 2 relevant for this study (the smallest 
scales probed by the kinematic data are R ~ 1.5 kpc ~ rs/100), 
this parametrization accurately describes the enclosed mass M{r) 
with an rms error of 0.1% and a maximum error of 1.8%, which is 
smaller than the uncertainties in cr^. For the Sersic profile, we fit 
16 Gaussian components to the surface density E(/2/ile) for many 
values of n, and then fit the amplitudes of the components with 4th 
degree polynomials in n. The result accurately describes Y1(R/Re) 
over the range 10“^ < R/Re < 5 and 3 < n < 7 relevant for this 
paper with an rms error of 0.3% and a maximum error of 0.9%. 


where imin is set such that q' > 0.5 as motivated by the 
absence of massive, non-rotating e llipticals with flatter 
shapes (e.g., Weijmans et al. 2014). We further assume 
that the DM halo follows the ellipticity and orientation 
of the stellar ellipsoids. This is motived by the observed 
close correspondence between the projected axis ratio 
and positio n an gle of mass and light in our lens mod¬ 
els (Section [T^ . 

Two free parameters other than inclination describe 
the stellar mass distribution. The first is an overall scal¬ 


ing 


asPS = Ty/Ty 


SPS 


(3) 


of the stellar mass relative to the SPS esti mate, which 


is cal led the IMF mismatch parameter (e.g., Treu et al. 
2010). As in previous work, this parametrization allows 


us to constrain systematic offsets in the stellar mass scale 
of the SPS models, which are premised on a Salpeter 
IMF. The second parameter is the anisotropy fiz of the 
velocity dispersion tensor in the meridional plane, for 
which we take a Gaussian prior iV(0.1,0.2) b ased on the 
sample of “slow rotator” ETGs considered in |Gappellari| 
er^((2007|. 

In addition to the asymptotic inner density slope /3, 
two scales are needed to describe the gNFW DM halo. 
Although the density profile is most easily written in 
terms of po and Ts, for comparison with simulations it is 
more useful to adopt the parameters M 200 and c_ 2 . Al¬ 
though the DM halo is ellipsoidal in our models, we define 
M 200 in terms of a spherical overdensity, i.e., the mass 
within the sphere of radius r 2 oo that has a mean density 
equal to 200 times the critical density. The concentra¬ 
tion c _2 = r 2 oo /?’-2 is defined in terms of the radius 
r _2 = (2 — /3)rs at which the local logarithmic density 
slope is —2, which is identical to the scale radius for 
NFW models. In order to properly compare with r 2 oo, 
we also use spherical measures of and r_ 2 jE 


Specifically, the sphericalized radii are defined as the geomet- 
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TABLE 7 

Mass Model Parameters and Priors 


TABLE 8 

Mass Model Eit Quality 


Parameter 

Description 

Prior 

log asps 

BGG Stars 

IMF mismatch parameter 

;7(-0.4,0.4) 

VTy 

dlogTfPS/j^iog^ 

iV(-0.15,0.03) 

Pz 

JAM anisotropy 

Af(0.1,0.2)_ 

COS i 

inclination 

F(0, imin) (Sjejl 

p 

DM Halo 
inner slope 

17(0,2) 

log AL200 

halo mass 

Table 3 

logC_2 

concentration 

17(0,1.3) 

^ext 

Other 

external convergence 

See iJO 
Af(1.0,TrU5) 

ffVD 

cr calibration 


Note. — U{a, b) denotes a uniform prior over the interval [a, b], 
and a) denotes a Gaussian prior with mean /i and dispersion 
a. A Gaussian prior is placed on logMgn^o based on the measure¬ 
ment from satellite kinematics in Table 1^ for the two groups that 
lack such a measurement, we take [7(133^, 14.5). 


For a given set of model parameters, we use JAM to 
construct the line-of-sight second velocity moment 
on a grid of radii and azimuthal angles. The surface 
brightness / and /crfog are then interpolated onto a rect¬ 
angular grid, smoothed by the PSF of the spectroscopic 
observations, and binned in the same rectangular aper¬ 
tures used for spectral extraction (accounting for the slit 
width and orientation listed in Table [^. We thus gener¬ 
ate a set of model velocity dispersions <7™°'^, which are 
compared with the data 0 '°*’®: 


XYD=y] 


A? 


(4) 


where is the uncertainty in measurement i. gvD is a 
calibration factor that accounts for correlated systematic 
uncertainties in the velo city dispersion measures. Based 
on the tests in Section |3.4[ we place a Gaussian prior 
A^(1.0,0.05) on gyu- 

The mass within a cylinder of radius 0Ein is then com¬ 
puted and normalized by the critical surface density for 
lensing to obtain the model convergence K^oup • This is 
compared to the measurements in Table 

) 2 

(^) 



Here is measured from the lens models with un¬ 

certainty (T„, and Next is the external convergence arising 
form foreground or background structures, which is de¬ 
scribed separately below. 

The model likelihood is then L oc exp[—(xy^ -|-x^)/2]. 
We incorporate the measurement of i^oo from satellite 
kinematics, where available (see Table M , via a Gaussian 
prior. The parameter space is explorecTusing the MGMG 
engine MultiNest. Table summarizes the parameters 
of our model and the priors. 


6 .1. External Convergence 

Several lenses show signs of foreground or background 
structures, either through secondary peaks in the redshift 


ric mean of the axes of the isodensity ellipsoid. 


Name 

XvD /-^VD 

xl 

Xm200 

GSWA107 

25.76/10 

0.00 

1.11 

GSWA141 

17.07/9 

0.85 

1.06 

GSWA164 

3.05/5 

0.47 


GSWA165 

7.61/5 

0.05 

0.08 

GSWA6 

5.48/6 

0.04 

0.05 

GSWA7 

2.73/7 

1.86 

1.05 

EOCL 

6.85/7 

0.18 

0.64 

J09413 

5.98/7 

2.06 

0.47 

GSWA163 

0.92/3 

1.46 

0.01 

GSWAl 

2.83/7 

0.86 



Note. — metrics are shown for the maximum a pos¬ 
teriori probability model found using our fiducial modeling 
procedure. A^vd is the number of bins in the ct{R) profile. 


distribution (Section |3.2[ ) or through external she ar Text 
required by the strong lens models (Section |4.2[ ). The 
three systems with the highest Text (GSWA7, J09413, 
EOCL) correspondingly show clear secondary peaks in 
the redshift distribution. GSWA6, which is ht with low 
Text) does have secondary redshift peaks, but they are 
located far (> 1 Mpc) from the lens. Reassuringly, of 
the three systems where the image configuration or data 
quality precluded a reliable estimate of Text (GSWA107, 
CSWA141, GSWA163), none shows a second redshift 
peak. Thus, the redshift survey and strong lensing anal¬ 
ysis are broadly consistent indicators of the presence of 
signihcant external structures. 

Strong lensing measures the mass within the Einstein 
radius, including any contribution from external conver¬ 
gence Next- To account for this, we introduce Next as an 
additional parameter. Eor the 6 systems with a small ex¬ 
ternal shear (Text < 0.04) or no measurement, we place 
a prior on Next of A(0, 0.04) based on the sm all fluctua¬ 
tions expec ted from large-scale structure (e.g., Takahaslri] 


et al.|2011 ). For the 4 systems with larger Text (GSWA7, 
JU9413, EOCL, CSWA165), we use Text to inform a prior 
on Next- If the external contribution is dominated by a 
single halo, the relation between Fgxt and Next depends 
on the slope of the density profile, with Text = «^ext in 
the case of an isothermal slope p oc We therefore 
take a log-normal prior on Next with a mean of Text but 
allow a broad dispersion of a factor of 2. An alternative 
approach is to estimate the distribution of Next condi¬ 
tional on another lens observable, such as Text or the 
local galaxy overdensity, using cosmological simulations. 
For galaxies with Text ^ 0.1, typical of our sample, this 
method yields yields a similar or slightly smaller disper- 


sion in , giving 

confidence that our prior is reasonable 

(e.g., 

Oguri et al. 

2005 

Suyu et al.||2010|. 


6 .2. Fit Quality 

Our mass models are generally flexible enough to fit all 
of the available data acceptably, as Table [^demonstrates 
via metrics. The main deficiency is that the steep 
rise in the velocity dispersion profiles of CSWA141 and 
CSWA6 is not well fit (Figure j^. This might signal that 
our mass models are not fully adequate in these cases. 
For the remainder of the paper, we bear this in mind 
when interpreting our results and pay attention to the 
influence of these two systems. 
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6.3. Consistency Between Mass Models and Lensing 

We have chosen to include only 0Ein as a lensing-based 
constraint in our mass models. However, strong lens¬ 
ing carries more information about the mass distribu¬ 
tion. In particular, the radial magnification can con¬ 
strain the logarithmic slope of the density profile 7 sl (see 
Table [^. This slope is affected by the well-known mass- 
sheet aegeneracy and is therefore less robust than the 
mass within 0Ein- Nevertheless, it is interesting to test 
whether our mass models are consistent with the avail¬ 
able strong lensing data in a broader sense. With this 
aim we compared 7 sl ) which was derived solely from the 
strong lensing data (Table [^, to the total density slope 
7 LD,eBin of ffio mass models introduced in this section, 
which are fit to the stellar kinematics and 0Ein- Since 
strong lensing is sensitive to the mass profile over the ra¬ 
dial range where multiple images are formed, we measure 
7LD,eEi„ locally around R = 0Ein- CSWA7 was excluded 
from this comparison since the pixel-level data were not 
used in its strong lensing fit (see Section 6). 

The mean density slopes recovered from the two meth¬ 
ods agree very well: (7sl) = 1.59±0.08 and (7LD,eEi„) = 
1.57±0.05. On an object-by-object basis, the scatter be¬ 
tween the two methods is = 0.28. Accounting for the 
statistical errors in 7LD,eEin estimated from the Markov 
chains, statistical agreement between the methods would 
require errors of ~ 0.2 in ygL. Although it is difficult to 
independently evaluate the uncertainties in 7sl, which 
are dominated by systematics, we note that t he value of 
7 ST, de rived for CSWAl in a detailed study by |Dye et al. 
(2008) is 0.30 higher than ours, possibly due to differ- 


^- 0.47 


ences in the source model and the subtraction of the lens 
galaxy and background. This comparison demonstrates 
that the mass models developed in this section have den¬ 
sity slopes consistent with the strong lensing data for 
reasonable estimates of the uncertainties in both meth¬ 
ods. 

The largest outlier in this comparison is CSWA107, 
with 7LD.eEi„ = 1-83 ± 0.19 but 7 sl = 1.18. Since a 
shallow projected density profile is preferred to produce 
the naked cusp configuration, the disagreement in this 
system may reflect a degeneracy with the ellipticity, ex¬ 
ternal shear, or the effect of the mass-sheet degeneracy. 
The good agreement between 7 sl and 7LD,eEi„ on aver¬ 
age, however, shows that the mass-sheet degeneracy is 
generally not a limiting factor on these scales. 

7. THE MASS-CONCENTRATION RELATION AT THE 
GROUP SCALE 

The various data sets assembled in this paper con¬ 
strain the group mass distribution on multiple scales. In 
this section, we consider the large-scale mass distribution 
as quantified by the halo mass-concentration relation 
(MCR) and compare it to theoretical relations obtained 
from CDM simulations. Here we consider only the 8 
lenses in our sample for which we have measu red t he halo 
mass M 200 from satellite dynamics (Section 3.3). Since 
theoretical studies of the MCR are usually premised on 
an NFW DM profile, our results in this section only 
are derived using mass models that assume NFW ha¬ 
los (i.e., /? = 1 and c _2 = C 2 oo)- Furthermore, be¬ 
cause the MCR evolves with redshift, we evolve all mea¬ 
sured concentrations—including those drawn from the 
literature—to the mean redshift zq = 0.36 of our sample 


using the Duffy et al. ( 2008|_^caling C 200 oc (1 ± z 

The left panel of Figure]^ shows our constraints on 
the halo mass and concentration of each lens. Lines of 
constant 0Ein are overlaid. As expected, the covariance 
between M 200 and C 200 follows the slope of these lines, 
since 0Ein is more precisely measured than M 2 oo- Con¬ 
straints are broad for individual lenses, primarily due to 
the uncertainties in the satellites’ velocity dispersion and 
hence M 2 oo- Therefore, we combine results from the 8 
groups using a hierarchical Bayesian method that allows 
us to infer the mean (logM 2 oo), the mean concentration 
(logC 2 oo) -^200 = 10^'* M©, and the intrinsic scatter 
in both distributions. The mathematical details of this 
framework are given in Appendix B. 

The resulting constraints are (logM 2 oo) = 14.0 ± 0.1 
and (logC 2 oo) = 0.70±0.07 (thick black ellipse in the left 
panel of Figure [9). This mean concentration is entirely 
consistent with the theoretic al expectation for u nmod - 
ified NFW halos derived by Dutton & Maccio (2014), 
who predict a mean log C 200 = 0.67 at log M 200 = 14 and 
z = 0.36. This concentration is slightly higher than some 
earlier theoretical MCRs, wh ich[Dutton fc M accio (|2014|) 
attribute to their use of the FlanctTColiaboration et al.l 


(2014) cosmological parameters. To illustrate the differ¬ 


ence comp ared a theoret i cal M CR based on the WMAP5 
cosmology, |Duffy et al.l I 
ds 


(2008) found logC 2 oo = 0.56 at 
lift when considering all halos 


the same rnass and reds 
and logC 2 oo = 0.62 when restricting to the relaxed ones. 
We also find that the intrinsic dispersion in concentra¬ 
tions, CTiogc' = 0.14 ± 0.07, is consistent with the range 
predicted in CDM simulations. 

Lens selection can have an important effect on the con¬ 
centrations of a sample. The selection function of surveys 
that search for strong lenses is complex, but the Ein¬ 
stein radius 6*Ein is the s ingle most important variable 
(e.g., Gavazzi et al.||2014 ). As Figure j^shows, contours 
of constant ^Ein are diagonal in the mass-concentration 
plane. This implies that selecting lenses within a par¬ 
ticular range of 0Ein (or with a non-uniform weighting) 
will lead to biases in the slope and intercept of the MCR 
relative to that of the underlying halo population, un¬ 
less the selection is taken into account. We estimate this 
effect using a halo o ccup ation model. The details are in¬ 
troduced in Section [ 9 ( 1 ] but for the present purpose it is 
sufficient to consider ETCs resi ding in a cosmological dis- 
tribution of halos following the Dutton & Maccio (2014) 
MCR. Erom such a mock sample, we weight galaxies to 
match the distribution of 0Ein in our group lens sample 
and compute the mean halo mass and offset from the 
MCR. We find (log M 200 ) = 14.0, in agreement with our 
dynamical measurement, and (AlogC 2 oo) = O.llj^ In 
other words, lenses with 275 < 0Ein < 571 follow a^CR 
that is somewhat offset from that of the parent ETC pop¬ 
ulation. Applying this estimated correction to our infer¬ 
ence for the group lenses yields (logc^oo) = 0.59 ± 0.07 
for the mean concentration of the underlying halo pop¬ 
ulation at M 200 = 10^"^ Mq, which is shown by the thin 
ellipse in Figure [9] 

This correctionfor the Einstein radius selection has the 
same magnitude as the differences among current theo- 

Here we populate NFW halos with galaxies having a Salpeter 
IMF, b ut va riations to the IMF and inner DM profile, described in 
Section |9lT| affect this correction only at the ~ 0.03 dex level. 
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Fig. 9.— Left: Halo masses and concentrations normalized to zq = 0.36. Colored contours show the 68% credible regions for individual 
lenses. The thick black ellipse shows our hierarchical inference for the mean halo mass and concentration of the group lenses and their ztlcr 
uncertainty. The thin ellipse represents our estimate of the underlying halo population after accounting for our selection in Einstein radius. 
These are compared to t he theoretical relations indicated in the lower caption. Dotted black line s ind icate the ztlc scatter in logC2oo from 
[Dutton &: Maccio| l |2014l l. Contours of ^Ein derived from our halo occupation model (Section jQ.lj these are insensitive to details of the 
stellar distribution for ©Ein ^ 2".) Right: Comparisons to published MCR constraints around ~ 10^^ Mq. The individual jOguri et al.j 
| |2012[ l lensing mea surements a.nd the lBuote et al.lll2007ll X-ray data are binned to produce the blue and gray squares, respectively. Both 
the stack of the fulljFoex et al.j (|2014|| weak lensing sample (solid red circle) and their stacks in three bins of arc radius (red error bars) 
are shown: from low to liigli C 200 , 2" < R^ < 3('5, 3('5 < Ra < 5'.'5, and 5'.'5 < Ra < 20". 


retical MCRs (~ 0.1 dex). We conclude that, within the 
present uncertainties, our group-scale lenses have a mean 
concentration consistent with unmodified, cosmologically 
motivated halos. 

7.1. Comparison to Published Concentrations 

Since the concentrations of groups and low-mass clus¬ 
ters are contentious, here we briefly compare our results 
to other studies. We include only those that extend as 
low as M 200 — IO^'^Mq, convert published masses and 
concentrations to our overdensity definition, and evolve 
concentrations to zq = 0.36 using the scaling described 
above. Since selection effects will prove to be important, 
it is useful to have a model for the relative numbers 
of strong len ses having different Einstein radii. 


More 


et al. (2012) showed that the image separation distri¬ 
bution (ISD) roughly follows dP/dO oc 0“^ ® for image 
separations 9 ~ 20Ein = 3" — 30". 

Some authors studying group-scale strong lenses have 
found “nor mal” concentration s consistent with our mea¬ 
surements. [Auger et al.|(|2013|) analyzed 26 strong lenses 
(median ^Em = 4'.'0j selected from the CASSOWARY 
catalog, the source of many of the lenses in our group 
sample. They combined 0Ein with an estimate of M 200 
for each lens based on galaxy richness, using a scaling 
relation calibrated to X-ray massesj^ The right panel 
of Figure [^ shows that the Auger et al. concentration 
at M- 2 ()o = 10 ^ ^ Mfn (green star) agrees with our value. 
Deason et al.j (|2013|) analyzed CSWA 163 using similar 
methods to our own. As mentioned in Section we 
have incorporated their data into our sample, and we 
find consistent values of M 200 and C 2 oo- 

Comparing the richness-based masses with our dynamical 
masses for the 7 systems in common, we find agreement within the 
uncertainties in 5 cases. The richness-based masses are much higher 
(0.8 dex) for EOCL—likely explained by the presence of multiple 
structures at similar redshifts artificially boosting the richness— 
and CSWA141. 


Other authors have claimed evidence for overconcen¬ 
trated halos at the scale of groups and low-mass clusters. 
In an thorough analysis of 28 le nses discovered in the 
Sloan Giant Arcs Survey (SGAS), Oguri et al. (2012) in¬ 
ferred a very steep slope for the MGR, Cvir oc 
(c.f. Cvir oc in CDM simulations). Oguri et 

al. modeled two selection effects: the probability for a 
cluster to produce an arc with a length-to-width ratio 
l/w > 5 (i.e., the arc cross-section), and the selection 
function of the SGAS among suc h len ses, which they ap¬ 
proximated as proportional to V^Ein ■ Since these effects 
were insufficient to reconcile GDM-only simulations with 
the steep observed slope, Oguri et al. suggested that ha¬ 
los are significantly modified on large scales by baryon 
cooling, even at fairly high masses > IO^'^Mq. We sug¬ 
gest that the steep MGR found by Oguri et al. arises 
from the SGAS sample selection. In particular, there are 
almost no lenses with 0Ein < 5" (see blue points in Fig¬ 
ure right panel), even though the number of strong 
lenses in the universe incre ases rapidly toward smaller 
^Ein (e.g.,|More et al.||2012|. This implies a rather hard 
cutoff in tne selection function. As Figure shows, this 
corresponds to a diagonal cut in the mass-concentration 
plane that will induce a steep slope, which matches that 
of the Oguri et al. MGR almost exactly. Samples of lenses 
with ^Ein > 5" are adequate to measure the concentra¬ 
tions of massive clusters with M 200 ^ Mq, but one 
must probe smaller 0Ein to reach halos with “normal” 
concentrations and smaller masses. Similar considera¬ 


tions apply to the Wiesner et al. (20121 analysis of the 
S loan Bright Arcs S urvey sample. 


Foex et al. (2014) used a stacked weak lensing anal- 
ysis to measure the mean mass and concentration of 
80 strong lensing groups in the SARGS sample having 
2" < 0Ein ^ 20" (mean ~ 4"). Their mean halo mass is 
comparable to our group sample, but their mean concen¬ 
tration is significantly higher (orange circle in Figure 
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TABLE 9 

Total Density Profile Slopes and Stellar Fractions of 
ETG Lenses in Massive Halos 


Name 

Re (kpc) 

logM*/M0 

7tot 

f2D 

*,Salp 


Group-scale Lenses (BGGs) 

CSWA107 

25.8 

11.95 

1.63 ±0.19 

0.16 ±0.07 

CSWA141 

40.3 

12.03 

1.50 ±0.11 

0.09 ±0.03 

CSWA164 

10.8 

11.85 

1.84 ±0.13 

0.30 ± 0.10 

CSWA165 

29.6 

12.10 

1.62 ±0.10 

0.15 ±0.05 

CSWA6 

30.0 

12.13 

1.52 ±0.12 

0.16 ±0.05 

CSWA7 

24.2 

11.96 

1.62 ±0.15 

0.18 ± 0.05 

EOCL 

25.7 

12.03 

1.68 ±0.14 

0.19 ± 0.05 

J09413 

46.1 

12.26 

1.57 ± 0.15 

0.12 ±0.03 

CSWA163 

17.6 

11.89 

1.74 ±0.13 

0.19 ±0.03 

CSWAl 

20.9 

12.05 

1.66 ±0.10 

0.16 ±0.02 

Clusters-scale Lenses (BCGs) 

A611 

52.0 

12.43 

1.26 ±0.04 

0.06 ±0.01 

A383 

44.0 

12.30 

1.14 ±0.08 

0.06 ±0.01 

A2667 

54.0 

12.15 

0.96 ±0.08 

0.04 ± 0.01 

MS2137 

31.0 

12.23 

1.25 ±0.05 

0.07 ± 0.01 

A963 

34.0 

12.32 

1.27 ± 0.07 

0.10 ± 0.02 

A2390 

28.0 

12.06 

1.34 ±0.05 

0.07 ± 0.01 

A2537 

59.0 

12.48 

1.02 ±0.04 

0.05 ±0.01 


Note. — All stellar masses refer to a Salpeter IMF. For the 
BGGs, Rg is measured in rest-H band and M, = Ly (Ta¬ 

ble |m. The BCGs have been fit with de Vaucouleur’s profiles and 
so oirFer from published and M* in |Newman et al.l (|2013^. 
Random errors in Mt/Lgps are estimated as U.U5 dex (Section [y 
for the BGGs and 0.07 dex (N13a) for the BCGs. The random 
error in includes both the projected mass and M*/LsPS- 


right panel). By stacking in bins of arc radius Ra, a 
proxy for they note a strong increase in concen¬ 

tration, but not halo mass, with increasing Ra (orange 
crosses in figure). The concentration in their lowest bin 
2 " < Ra < 3'.'5 is consistent our measurement, suggest¬ 
ing that the cause of the discrepancy in the full stack 
can be traced to the lens sample, not a difference in the 
mass probes used. Although the origin of the difference 
in concentration relative to our sample is not fully clear, 
it is possible that higher-0Ein lenses are somewhat over¬ 
represented in the Foex et al. stacks 

Finally, it is interesting to compare to results from X- 
ray studies. Figure shows that for halo masses M 200 ^ 
10 ^^ Mq, the Buote et al. ( 2007 1 compilation agrees with 
both our group sample and the Dutton & Maccio (20141 
theoretical MCR. The Buote et al. concentrations exceed 
the theoretical relation only at lower masses. 

In summary, we find “normal” concentrations in 
10 ^^ Mq halos indicating that baryons have little effect 
on the DM distribution around the scale radius. This 
echoes some earlier work at this mass scale, but not all. 
Differences in lens sample selections are li kely to explain 
at least some of the differences. In Section [TO] we discuss 
our results in the context of hydrodynamical simulations. 


8 . THE TOTAL MASS DENSITY PROFILE WITHIN THE 
EFFECTIVE RADIUS 


Approximately equal numbers of lenses with Ra = 2" — 375 
and Ra = 375—575 are present in their stacks, whereas the More et 
al. ISD implies that the former should be 3 times more numerous. 


Moving to smaller scales, we now consider the den¬ 
sity structure within the effective radius. We begin with 
the average logarithmic slope ytot of the total density 
profile within Re, where ptot oc This is one of 

the simplest and most robust quantities that can be in¬ 
ferred from a leasing and dynamics analysis. Further¬ 
more, since it requires only a measure of 0 Ein and a sin¬ 
gle aperture velocity dispersion, ytot has been measured 
for large samples of galaxy-scale lenses and so is partic¬ 
ularly useful for examining trends across a wide range of 
ETG properties. In this section we combine our group- 
scale lenses with earlier data on galaxy- and cluster-scale 
lenses. We examine empirical trends in ytot within and 
among these samples, which collectively span a factor of 
~ 60 in halo mass. In Section we will interpret these 
trends using CDM-motivated models and focus on the 
more challenging goal of separating the luminous and 
DM density profiles. 

Since our mass models have separate DM and stellar 
components, the total density pr ofile is not explic i tly pa - 
rameterized. Instead, following Dutton & Treu (2014), 
we define ytot as the mass-weighted mean density slope 
within Re'. 


iTTRlp{Re) 
M{Re) ’ 
( 6 ) 

where p(r) and M{r) are the spherically averaged density 
and enclosed mass profiles. The posterior distribution of 
7 tot for each lens can then be evaluated from its Markov 
chains. Heuristically, the data constrain the density slope 
most robustly between the Einstein radius and the ve¬ 
locity dispersion aperture. Since the Einstein radius is 
typically ~ 0.5 — li?e (see Figure [^, Re is a convenient 
location for defining ytotP^ 

To compare our group-scale lenses with other samples 
of ETG lenses, we have gathered data from the following 
sources: 




Ttot — 


M{r) 


A-kt^ pir) 


dlogp 

d\ogr 


dr = 3 


• Galaxy-scale lenses: We select 59 lenses from the 
SLAGS survey wit h measured Rs, stellar masses, 
and density slopes (Auger et al. 2009 2010a). Since 
Auger et al. use single-component power law mass 
models, their slope 7 ' is nearly equivalent to our 
7 tot (see Sonnenfeld et al.|[2015 ). The average 0Ein 
is 174 (rescaled here, as throughout, to a leasing 
distance ratio D^s/Ds = 0.7). Although the halo 
mass cannot be measured for individual SLAGS 
lenses, a stacked weak leasing analysis of a sub- 
set indicated a mea n {\ogM 2 oo/MQ) = 13.2 ± 0.2 
(IGavazzi et al.|2 007l converted to ou r adopted over- 
density; see also ^ger et al.||2010b ). Based on the 
distribution of 0 Ein and our halo occupation model, 
introduced below, we estimate the scatter in M 200 
within this sample to be ^ 0.3 dex. 


• Group-scale lenses: We use the 10 lenses analyzed 
in the present paper having 0Ein = 275 — 571. 
As shown in Section the average halo mass is 

While this represents a slight extrapolation for our group sam- 
pie, where ^Ein — 0.7-Re on average, quantitative differences in the 
slope defined over slightly different radial ranges are small. For 
example, the mean (7tot} would change by only 0.06 if we were 
instead to define it within 0.7Re- 
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Star formation efficiency 



6»Ein [arcsec] log s, [Mg kpc^^] 


Fig. 10.— Trends in the total density slope 7tot and projected stellar mass fraction within Re are plotted as a function of Einstein 

radius ^Ein and mean stellar surface density S* = M*/{27rR‘l) for three ETG lens samples spanning 10^^ Mq — 10^^ Mq in halo mass: 
galaxy-scale lenses drawn from the SLAGS survey (blue circles), the present group-scale sample (black squares), and the BCGs of massive 
clusters studied by N13 (green stars). Darker symbols show the mean of each data set. Left panels: Sold lines indicate the mean trends 
for the 5 halo occupation models introduced in Section |9.1| which adopt different IMFs and inner DM density profiles. For comparison, 
the dashed “NFW only” line corresponds to an unmodified NFW halo with no stars. The top axis delineates the average star formation 
efficiency M*/M 2 oo/(^£i/^m) as a function of ^Ein; which is taken from the Salpeter-|-NFW halo occupation models. Right panels: In 
three bins of halo mass, lines indicate the trend in the Salpeter-|-NFW halo occupation model along with its Icr scatter (dashed). 


{log M 200 /Mq) = 14.0 ± 0.1 with an intrinsic scat¬ 
ter of 0.2 dex. The relevant measurements for the 
group- and cluster-scale lenses are listed in Table 


• Cluster-scale lenses: For the central galaxies of 
massive clusters , we use t he 7 BCGs analyzed by 
Newman et al. (2013a|b). In order to compare 


more consistently with the other lens samples, we 
remeasured the BCGs’ Re and M* using a de Vau- 
couleurs’ profile (rather than the dPIE profile used 
by N13) and also recomputed ytot following Equa¬ 
tion]^ This yields a mean (ytot) = l.lSiO.OTto^o?! 
although this definition of 7tot differs formally from 
that adopted in N13a, in practice the difference is 
quite small (c.f. (qtot) = 1-16 in N13a). The BCGs 
have a typical 0Ein = U" and occupy halos with a 
mean mass of (logM^m/^Q) = 14-9 and a scatter 
of 0.3 dex. 


The group sample is situated at only slightly higher 
mean redshift {{z) = 0.36) than the SLAGS {{z) = 0.20) 


and N13 ((z) = 0.25) samples, so evolutionary differences 
among the samples are expected to be minimal. 

In addition to ytot, we measure the projected stellar 
fraction within Rp for each lens: 


/*,Salp = Mt^{Re)/M,i^,iR, 


(7) 


where M^^{Re) is the projected mass within Rp derived 
from the lensing and dynamics model, and M^^{Rp) is 
the stellar mass in the same aperture estimated using 
SPS models and a Salpeter IMF. 


Figure 10 illustrates the trends in qtot and /* s^ip 


among the three observational samples as a function 
of 0Ein (left panels) and stellar surface density S* = 
M*/(27ri?g) (right panels). The first trend clearly visible 
in the data is the decline in ytot and /* with increas¬ 
ing 0Ein and M 200 - The mean total density slope in the 
group-scale lenses is (ytot) = 1.64±0.05±0.07 (the second 
error is our estimate of the systematic uncertainty; see 
Section 9.3). This is significantly shallower than that of 
the galaxy-scale lenses, (qtot) = 2.09 ± 0.03, and steeper 
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than that of the BCGs, (7tot) = 1-18 ± O.OT^Q gy. The 
intrinsic scatter within all samples is much smaller than 
the systematic variation among them: a-y = 0.07 ± 0.05 
(groups), 0.19±0.03 (galaxies), and 0.17±0.08 (clusters). 
As we will see, these trends arise from a dependence of 
7tot and /* g^jp on halo mass. The relatively small scat¬ 
ter within each sample arises from the similarly narrow 
ranges of halo mass that they span, while the system¬ 
atic trend among samples arises from the fact that they 
occupy nearly disjoint ranges of 0Ein and M 2 oo- 

The stellar fraction also declines sharply with increas¬ 
ing Einstein radius and halo mass: (/* gaip) = 0-60 for 
the galaxy-scale lenses, 0.17 for groups, and only 0.06 for 
clusters. To first order, this trend drives the decline in 
7tot: since the DM distribution is more extended than 
the stars, 7tot will be shallower in systems with lower 
/* Saip- The main question is whether the trend in 7tot 
can be ascribed entirely to the declining stellar fraction, 
or whether a non-universal IMF or dark matter profile is 
required. We return to this question in Section 

In addition to the broad trends in 7tot ana /* g^ip 
evident from comparing the galaxy, group, and cluster 
lens samples with one another, there are also correla¬ 
tions within each sample. For the galaxy-scale lenses it 
has been found that E* = M^/{2ttR'1), the mean stellar 
surface density within i?e, is correlated with 7tot- This 
is natural, since galaxies with more concentrated stellar 
distributions should have steeper density profiles. Fur¬ 
thermore, the dependence on E* appears to be funda¬ 
mental, since th ere is no residual correlat ion with M* or 
Re individually ( Sonnenfeld et al.||2013b ). 

Our seco nd iniportant conclusion that we draw from 
Figure (right panels) is that although the central 
galaxies of groups and clusters have lower E* and shal¬ 
lower 7tot than the galaxy-scale lenses, they do not lie 
on a simple extension of the scaling relation seen in the 
galaxy-scale samples. The dashed line in the lower-right 
panel indicates the slope in ferred from the combine d 
SLAGS and SL2S samples by [Sonnenfeld et al.j (|2013b|); 
it clearly over-predicts 7tot for the group- and cluster- 
scale lenses. This implies that 7tot depends not only on 
properties of the stellar distribution, namely E*, but also 
depends explicitly on properties of the DM halo. 

To quantify this dependence and provide a reference for 
future observational studies and numerical simulations, 
we fit a linear regression 


7tot = 7tot,o + ^T 7 ^(log E* - 9.0) 


d log E* 
t^7tot 
d\ogM, 
97tot 


d log M 200 
d^tot 


(logM* - 12.0) 
(logM2oo - 14.0) 


+ -^{z-0.3) + N{0,a'^) 


( 8 ) 


to the combined sample of galaxy, group, and cluster- 
scale lenses. Here fV(0,cr(,) represents Gaussian intrin¬ 
sic scatter in 7tot after accounting for these linear de¬ 
pendences. A Bayesian procedure was used to infer the 
linear dependences of 7tot on each parameter. To ac¬ 
count for small redshift differences among the samples. 


we t ake a Gaussian prior ./ V(— 0.30, 0.25) on d^tot fdz 
from Sonnenfeld et al. (2014|), but otherwise we use uni¬ 
form prldrsl Individual halo masses are unavailable for 
the SLAGS lenses. Therefore we introduce a parame¬ 
ter (logM 200 , slags) with a Gaussian pri or A^(13.2,0.2), 
based on stacked weak leasing results (Gavazzi et al. 
20^ ; an uncertainty of 0.3 dex is assigned to the halo 
mass of each SLAGS lens based on our estimate of the 
intrinsic scatter]^ We find the following linear depen¬ 
dences for 7tot: 

70,tot = 1.79 ±0.06 

djtot 


d log E* 

^7tot 

91ogM* 

^7tot 

d log M 200 


= 0.34 ±0.11 
= -0.06 ±0.12 
= -0.33 ±0.07 


a ' = 0.12 ± 0.10 


(9) 


(The posterior distribution for d"ftot/dz is consistent 
with its prior, as expected, and is omitted here.) 

We recover the trend with E* seen in galaxy- scale 
lenses with a consis tent value of i97tot/i9 log E* (|Son- 


nenfeld et al. 2013b|), as well as the lack of any addi¬ 
tional dependence on the stellar mass of the galaxy (i.e., 
57tot/91ogM* is consistent with zero). However, we also 
find a significant additional dependence on halo mass: 
d"/tot/dlogM 2 oo = —0.33 ±0.07 differs significantly from 
zero. The strength of this dependence on halo mass is 
roughly equal to the previously known one on E*. Earlier 
work has shown sy stematic variations in the density pro- 
file o f ETGs (e.g., ^nd et al.||2008[ Humphrey fc Buote 


2010| Newman ^ al.|2013a |Sonnenteld et al.|2bl3b ' Tor- 


tora et al. 2014p , but to our knowledge this is the first 

measurement of the bivariate dependence on two param¬ 
eters, one connected to the stellar distribution and the 
other to the DM halo. This is significant for the interpre¬ 
tation of t he “b ulge-halo conspiracy” as we will discuss 
in Section 10. 2| These findings are a further demonstra¬ 
tion of the non-homology of massive ETGs and show that 
their internal density structure cannot be predicted from 
the stellar light alone. 

One caveat is that the masses and radii of massive 
galaxies are sensitive to the light profile and measure¬ 
ment technique, especially for BGGs. To test the possible 
effect on our resul ts, we compared the Ey of our BGGs to 
measurements by Kravtsov et al. (2014), who took par¬ 
ticular care to fit the large-radius light profile of several 
BGGs including all intracluster light. Gonsidering their 
4 galaxies in halos with M 500 > 5 x 10^"* Mq, we find 
possibly lower E* by ~ 0.3 dex compared to our BGGs. 
Shifting the E* of our BGGs downward by this amount 
changes d^tot/d\og M 200 by only ±0.08. This is compa¬ 
rable to the random uncertainty and demonstrates that 
our basic result does not depend on the exact definition 
of the masses and radii of large galaxies. 


9. 


SEPARATING DARK MATTER AND STARS WITHIN 
THE EFFECTIVE RADIUS 


Likewise, for the two group-scale lenses lacking satellite kine- 
matic data, we use logM200 = 14.0 di 0.2 based on the mean and 
intrinsic scatter seen in the other 8 group-scale lenses. 
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log M./Mg 


Fig. 11. — Inputs to the halo occupation model described in Section |9.1| Left: The relation between halo mass a.nd the stellar r nass o f 
the central galaxy for three lens samples: galaxy-scale SLAGS lenses, based on the stacked weak leasing result of|Gavazzi et al.| ||2007|l; 
the present group-scale lenses, based on satellite kinematics in individual lenses; and the N13 clusters, based on iridividual weak lensing 
measures. Several theoretical relations are ov e rlaid. All stellar masses have been converted to a Salpeter IMF. The black line shows the 
abundance matching curve of [Kravtsov et al.| (|2014| Appendix A) shifted by —0.15 dex in logM* to better match the lens samples; this 
is the relation used in our halo occupation models. Black dashed lines show the la scatter at fixed M2oo- Right: The stellar mass-radius 
relation for lensing ETGs. A quadratic fit smoothly connects the galaxy, group, and cluster scale lenses. 


Although the total density profile within can be 
characterized well in individual ETG lenses (Section!^, 
the ultimate goal of separating the luminous and dark 
matter profiles is more challenging. This separation is 
well motivated because it entails (1) the distribution 
of DM on sub-galactic scales, which is sensitive to the 
poorly understood interplay between baryons and DM 
during galaxy formation (and possibly DM particle mi¬ 
crophysics; see introduction), as well as (2) an absolute 
mass scale for the stellar population, which is sensitive 


namics I 

Auger et al. 2010b| 

Treu et al. 

20101 Spiniello 

et al.|20] 

1.1[ Gappellari et al.j2 

013 

Newman et al.|2013b|) 


tion lines in the integrated galaxy light (jConroy fc van 


Dokkumi 2012| [Ferreras et al. 2013[ La Barbera et aT 


20Tm 


sive 


if Spiniello et al.|2014|) have both indicated that mas- 
F'i'Gs have a “heavier” IMF than the Milky Way, al¬ 


though the genuine degr ee of converg ence between these 
techniques is debatable (Smith|[2014). 


We take two complementary approaches toward dis- 
entangling the stellar and DM profiles. In Section |9.1 


we construct a simple, GDM-motivated halo occupation 
model that we use to forward model the distribution of 
7tot and /,^saip within the ETG population. By varying 
the model ingredients and comparing with the observa¬ 
tions presented in Section]^ we can constrain the small- 
scale DM distribution and the absolute stellar mass scale. 
This approach naturally encompasses the full range of 
galaxy- to cluster-scale ETGs and allows us to explore 
the origins of the trends discussed in Section]^ A disad¬ 
vantage is that it does not incorporate all of the data, in 
particular the resolved st ellar kinematics for our group- 
scale ETGs. In Section |9.2| we therefore compare the 
results from the first approach with those obtained from 
direct mass modeling using the full data set collected for 
our group-scale lens sample. 


9.1. Interpreting Density Profile Trends with Halo 
Occupation Models 

Here we introduce a set of models in which GDM- 
motivated halos are populated with ETGs that have, by 
construction, the same stellar masses and sizes as the 
ETG lens samples. Since the galaxy sizes and relative 
stellar masses are known empirically, they remained fixed 
throughout our set of models, whereas the inner DM pro¬ 
file and the absolute stellar masses (determined by the 
IMF) are varied. These ingredients are then constrained 
by studying the effect of their variation on the observ¬ 
ables 7tot and /*°aip- 

In shor t, we randomly samp le halos from a halo mass 
function (Angulo et al. 20121 at z = 0.36 and assume 
that they folfow a theoretical MGR. Galaxies are then 
assigned to the halos based on an empirical stellar mass- 
halo mass relation. Their stellar distribution is assumed 
to follow the de Vaucouleurs’ profile with R e sampled 
from the ob served stellar mass-radius relation. IDutton ^ 
Treu (2014) took a similar approach; the main difference 
with the present paper is that we consider trends over a 
wide range of halo masses, and we use scaling relations 
that are constructed from the lens samples themselves. 
Below we describe the main ingredients in more detail: 


Mass- concentration: We use the IDutton fc Macciol 
(2014) MGR at z = 0.36, the mean redshift of the 


group lenses, with a scatter of 0.13 dex in logC2oo- 

• Stellar mass-halo mass: The stellar masses and 
sizes of the model galaxies must be compatible 
with the lens samples. Since the luminosity and 
Re of massive ellipticals are sensitive to details of 
the measurement technique, we construct a stellar 
mass-halo mass relation using observations of the 
same lens ga laxies under study. In the left panel 
of Figure we compare t his to the relation pre - 
sented in tEe Appendix of Kravtsov et al. (2014), 
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Fig. 12.— Intrinsic scatter in 7tot f or t he same data sets and halo 
occupation models plotted in Figure with matching colors and 
symbols. For the halo occupation models, we select galaxies having 
a ^Ein within ±0.15 dex of the indicated value, which approximates 
the width of the ^Ein distribution in the lens samples. _[Note that 
cr-y is the total scatter and differs from in Equation!^) 

which w as constructed b y abun dance matching us¬ 
ing the Bernardi et al. (2013) stellar mass func¬ 
tion. These authors paid particular attention to 
photometry of massive galaxies. A small shift of 
—0.15 dex in logM, (after converting to a Salpeter 
IMF) brings this relation into reasonable agreement 
with the lens data. As emphasized by Kravtsov et 
al., many popular stellar mass-halo mass relations 
vastly underestimate M* at the high-mass end. We 
adopt a scatter of 0.2 dex in M* at fixed M2oo- 


• Stellar mass-radius: The right panel of Figure 11 
shows the relation defined by the lens samples. 
Since there is curvature in the relation at high 
masses, as noticed by earlier authors (e.g., Hyde 
& Bernardi 2009), we fit a quadratic relation: 
logi?e/kpc = 0.45 -I- 0.561ogMii -I- 0.26(logMii)^, 
where Mu = M*/10^^ Mq and the scatter in is 
0.11 dex at fixed M* □ 


We vary the stellar IMF between those of Chabrier 
(|2003|) and |Salpe ter| (|1955l), as defined in the |Bruzuai 
& Chariot (20U3() modeis, and consider three forms tor 
the inner i)M profile: unmodified NFW halos, ad i abati- 
cally contracted halos following the Gnedin et al. (2004) 
formulation, and cored NFW halos with a core radius 


?’core = 0.01r2oo (see Equation 3 of Newman et al. 2013a). 

For simplicity, we assume spherical symmetry and so 
do not model any orientation effects. For each mock 
galaxy we can then compute ^Eim 7toti and /fsaip- 
When computing mean properties of model galaxies, we 
weight the galaxies by the product of the halo mass func¬ 
tion (since the model galaxies are sampled uniformly in 


T he group-scale lenses are slightly offset from this fit (Fig¬ 
ure [13 because they are fit with free Sersic profiles, whereas the 
earlier galaxy- and cluster-scale lenses were fit with n = 4 fixed. 
This affects the comparison between data a nd models mostly in 
Salp slight offset for groups in Figure |lo[ upper-left panel) 
but has a small effect on the main focus of our analysis, 7tot: using 
n = 6 in the halo occupation models and increasing Re to better 
match the group-scale lenses shifts their (7tot} by only 0.04. 


logM2oo) and The second factor approximately 

accounts for the strong lensing cross-section, which is 
closely correlated with 0Ein (e-g-, Meneghetti et al.|2011 


although this factor has no effect when considering trends 
as a function of 6*Ein as we typically do). 

Lines in the lower-left panel of Figure show that 
all models predict shallower density profiles within in¬ 
creasing 0Ein or M 200 - To first order this decline in ytot 
is driven by the sharp decrease in /* g^jp visible in the 
upper-left panel, which is common to all of the halo occu¬ 
pation models. The detailed behavior of ytotj however, is 
sufficiently different to discriminate amongst these mod¬ 
els. Considering the galaxy-scale lenses first, we see that 
the mean (jtot) is best reproduced by the models with 
a Salpeter IMF. Both Chabrier-based models give too 
shallow slopes. This is consistent with earlier studies 
of the SLAGS sample that assume either unmodified or 


contr acted NFW halos (Auger et al. 2010b Treu et al. 
2010). The upper-left panel shows that the /* gaip data 


disfavor contracted halos with a Salpeter IMF, on aver¬ 
age, since the stellar fractions at small 0Ein are too low. 
Figure [^compares the scatter in qtot in the halo occupa¬ 
tion models to the intrinsic scatter measured in the three 
lens samples. The contracted Salpeter model predicts a 
scatter in ytot that is too small for the galaxy-scale lenses 
(although this scatter might increase if a distribution of 
parameters th at describe the adi abatic contraction were 
included, e.g., Gnedin et al. ]mT| ). Of the 5 halo occupa¬ 


tion models, therefore, those having a Salpeter IMF and 
unmodified or cored NFW halos best reproduce observa¬ 
tions of the ga laxy-scale lenses . This is consistent with 
recent work by Dutton & Treu (2014). 

Moving on to the groups, the mean qtot is consistent 
with either unmodified NFW halos and a Salpeter IMF, 
or contracted halos and a Chabrier IMF. This degen¬ 
eracy can be broken if we assume that the BGGs have 
a Salpeter-type IMF similar to the galaxy-scale lenses, 
which is reasonable if IMF variations are related to the 


central velocity dispersion (e.g., 

Treu et al. 

2010 

Cap- 

pellari et al. 

2013 

Posacki et al. 

2015) 

In that case. 


we would conclude that nearly unmodified NFW halos 
are favored, on average, since the slope of the contracted 
model is too steep, while that of the cored model is too 
shallow. Formally (ytot) is inconsistent with the con¬ 
tracted Salpeter model by 2.6tT and the cored model by 
4.2cr. We note that the /* g^jp data (upper-left panel) 
are not helpful for distinguishing models at the group 
and cluster scales. 

At the cluster scale, models with NFW or contracted 
halos have steeper ytot than all of the observations. This 
is consistent with N13, who inferred that the DM den¬ 
sity profile is shallower than the NFW profile within 
~ Re, whereas the total density slope 7tot is very close 
to that of a pure NFW halo (compare the green star 
and dashed line in Figure 10). Our model with a core 
radius of rcore = 0.01r2oo nts the clusters, which is un¬ 
surprising since it was intended to approximate the N13 
result. However, we see here that the same model does 
not match the group-scale lenses. This suggests that ei- 


If anything, the 0.06 dex higher average central dispersions 
of the BGGs relative to the SLAGS lenses would suggest a heav¬ 
ier IMF, although by A log M*/L < 0.1, which would place more 
tension on the adiabatically contracted models. 
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ther cores are confined to the most massive halos, or at 
least are smaller as a fractio n of r2oo in the groups. 

The right panels in Figure [l0| show that the halo mass 
dependence at fixed S* seen in the data is also present in 
the halo occupation models. Here we consider only the 
Salpeter+NFW model to clarify the trends, and we plot 
results (gray bands) in narrow intervals of halo mass that 
approximate those of the lens samples. At a given halo 
mass there is very little scatter in the relation between 
S* and 7tot or /* gaip) which again indicates that the 
range of ETG properties seen at fixed E* is driven by 
differences in halo mass. 

The halo occupation models allow us to estimate the 
effect that a strong lensing selection might have on the 
linear regression in Equation Although strong lens¬ 
ing ETGs are representative of ETGs having a given 
^Ein (or, nearly equivalently for the galaxy-scale lenses, 
a given cr), they are not strictly representative of ETGs 
of a given stellar or halo mass. We fit our ensemble of 
model galaxies using two different weights: (1) uniform 
in logM2oo, and (2) weighted to match the ^Ein distribu¬ 
tion of the full observed lens sample. Gompared to the 
uniformly weighted case, the lens analogs differ in 70, tot 
by -1-0.03, 97tot/i91ogE* by -1-0.07, and c^ytot/9 log M200 
by -1-0.03. These differences are less than the statisti¬ 
cal errors, demonstrating that for the range of galaxies 
we consider and the quality of the available data, the ob¬ 
served trends are minimally affected by th e strong lensing 
selection (see also Sonnenfeld et al.|[2015). 


9.2. Constraints on the Inner Dark Matter Profile from 
Resolved Stellar Kinematics 


T he f orward modeling approach presented in Sec¬ 


tion 


9.1 


relies on only two observables, ytot and /* g^jp, 
and so does not make full use of the data collected for 
our group-scale lenses, especially the resolved cr(i?) pro¬ 
files that are typically unavailable for galaxy-scale lenses. 
Here we consider the constraints on the mass distribu¬ 
tion of the group-scale lenses that arise from direct mass 
modeling (Section of the full data set presented in 
Sections 3-5. 

To quantify the inner DM profile, we introduce a mea¬ 
sure analogous to ytot- Specifically, we define 7 dm as 
in Equation but replace the total density and mass 
profiles p{r) and M{r) by those of the DM halo. This 
quantity represents the average DM density slope within 
and is better constrained by the data than the asymp¬ 
totic slope d that appears in the gNFW parametrization 
(Equatio n!^ . 

Figure [T^ shows the constraints on 7 dm and asPS for 
each lens. The lenses with the noisiest velocity dispersion 
profiles have broad contours indicating a nearly complete 
degeneracy between the DM profile and stellar mass nor¬ 
malization. As the quality of the data improve, they 
select preferred regions along this covariance. Broadly, 
the slope of the velocity dispersion profile is determined 
by the relative contribution of the stars and DM halo: 
the stellar mass profile contributes a roughly constant 
dispersion (p* oc r“^, a{r) oc r° for isotropic orbits) 
whereas the DM halo has a shallower density profile and 
so contributes a rising velocity dispersion at small radii. 
(Radially variable velocity anisotropy can also affect the 
shape of the velocity dispersion profile, and we explore 


this issue further below.) 

Since the decomposition between DM and stars is noisy 
in individual objects, we seek to combine the sample to 
constrain the mean agPS and 7 dm for the ensemble of 
group-scale lenses. To do so, we use the framework de¬ 
veloped in Appendix B and find (log agpg) = —0.11±0.06 
and (7 dm) = 1-35 ± 0.09 for our fiducial mass modelsj^ 
These constraints are displayed in Figure In ad(W 
tion to the constraints from our fiducial mass models, 
we also plot the systematic shifts that arise from several 
modifications of our model assumptions described fur¬ 
ther below. These are indicated by the colored symbols. 
The dashed blue contours indicate the posterior distri¬ 
bution if CSWA6 and CSWA141 are excluded, based on 
the poorer fits to their cr[R) profiles mentioned in Sec- 
tion |6.2[ To place these values in context, we select lenses 
with a matching dEin distribution in our halo occupation 
models and compute the mean 7 dm under the several 
models for the inner DM profile (dotted lines). We also 
include h ere a stronger adiabatic co ntraction prescription 
following 

Even a 


Blumenthal et al. (1986). 

:ter combining data from our 10 group lenses, 
there is a significant covariance between asPS and 7dm- 
This particularly true when considering the systematic 
shifts that dominate the error budget. Including these, 
we find (7dm) = 1.35±0.09 (stat.) I*!)];®! (sys.). Figure [l^ 
shows that asps is not robustly constrained by the group- 
scale lens data alone: it is particularly sensitive to the 
treatment of the radial M, /L gradient (compare the red 
and blue stars to the black circle) and to the inclusion of 
GSWA6 and CSWA141 (dashed contours). 

More informative results on the inner DM profile are 
possible with a prior on the IMF motivated by exter¬ 
nal observations. Several studies, i nclu ding our inference 
for group-scale lenses in Section |9.1| and for BCGs in 
Newman et al. (2013b), favor a Salpeter or even heav- 
ier llVIh' in massive ellipticals (see also references in Sec¬ 
tion]^. Therefore, we also consider the result of impos¬ 
ing a mea n S alpeter IMF in our analysis: the right panel 
of Figure 14 shows the constraints on (7 dm) adopting a 
Gaussian prior A^(0.0,0.05) on (logasps)- 
Under this assumption, the data favor an inner DM 
slope between that of an unmodified NFW hal o and a 


mildly cont racted one (i.e., contraction followin g Gnedin 
et al. 2004). Stronger Blumenthal et al. (1986) contrac- 


tion IS excluded, and the cored INh'W profile preferred 
by the clusters (rcore/T’200 = 0.01) is disfavored. How¬ 
ever, the significance of the latter conclusion is sensitive 
to the treatment of CSWA6 and GSWA141 (light blue 
line). When these systems are omitted, the discrepancy 
with the cored NFW profile decreases from 4.4cr to 1.9cr. 
These results are consistent with the conclusions drawn 


from our halo occupation models in Section 9.1 although 
the levels of significance differ. The present approach is 
more general, since it imposes no strong prior on halo 
concentrations and allows for galaxy-to-galaxy scatter in 
asPS- For this reason, we conclude that although the 
group-scale lenses disfavor the shallow DM profiles found 
in the centers of clusters, this difference is not yet defini¬ 
tive. 


Here osPS is defined relative to a Salpeter IMF and so is lower 
than a Chabrier-based definition, as used by N13, by 0.25 dex. 
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Fig. 13.— Credible regions for the IMF mismatch parameter asPS ^rid 7 DM 9 tti® average DM slope within Re 
lenses. Contours enclose 68% and 95% of the posterior distribution. 
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Fig. 14.— Left: Constraints on the mean (logasps) {7dm} for our ensemble of 10 group-scale lenses. The black circle and contours 
correspond to the fiducial mass model. Colored symbols show the systematic shifts that arise from various modifications to these model 
assumptions as indicated in the caption. Contours enclose the 68% and 95% credible regions. Right: Marginalized constraints on the mean 
DM slope (7 dm) of fh® group lenses if a prior centered on a Salpeter IMF is adopted for (log cusps). Lhe thick black line represents our 
fiducial mass model, while colored lines have the same meaning as the corresponding symbols in the left panel. Labels on the top axis, as 
well as dotted lines in the left panel, show the 7nA/r expected for various models of the inner DM distribution, derived by selecting lenses 
from our halo occupation models (Section |9.1[ ) matching the ^Ein distribution of the group lenses. We emphasize that 7 dm is averaged 
within Re and is not the same as the asymptotic inner slope /3 appearing in the gNFW parametrization (Equationly ; in particular, 7 dm 7 ^ 1 
for NFW halos. 

9.3. Systematic Uncertainties 

Here we outline several tests designed to assess the 
sensitivity of (7 dm) and (log asPs) to our model assump¬ 
tions. Colored lines and symbols in Figure show the 
inference obtained under each modification to the fidu¬ 
cial model described below. 

First, we switch from axisymmetric to spherical dy¬ 
namical calculations, as used in most dynamical analyses 
of lenses, and adopt a prior of iV(0.0,0.2) on the radial 
anisotropy paramet er f3r based on the results of |Cap-| 


derived by different a uthors (e.g., |Gerhard et al.||2001 
Cappellari et al.||2007|) , and a model with radiall y vary¬ 
ing anisotropy lollowing the Mamon & Lokas (20051 


pellari et al. (20071. Assuming that the axisymmetric 
models are rnore correct, the simplifying assumption of 
sphericity tends to make the IMF normalization slightly 
too light and the DM profile slightly too steep, but the 
effect is small (see Figure 14). Second, we considered 
spherical models with fixed velocity anisotropy: both ra¬ 
dially constant models with ,5,. = 1 — (T|/cr^ = ±0.25, 
which is intended to approximate the systematic differ¬ 
ences in the mean anisotropy of massive ellipticals as 


parametrization with = 0.5J?ek' The radially variable 
model gives results very close to the spherical isotropic 
model. These tests are reassuring for N13 study of BCGs 
based on spherical dynamics and constant anisotropy 
models. 

The biggest systematic uncertainties come not from 
the dynamical modeling, but from the magnitude of 
the M*/L gradient and the influence of CSWA6 and 
CSWA141. The red star in Figure show the results 
using our fiducial dynamical model but a radially invari¬ 
ant Ty, which we take as the mean (T|?®) within i?e 
listed in Table Neglecting the M*/L gradient entirely 
has a very large effect on (logasps), shifting it systemat- 


The anisotropy parameter increases from /Jr = 0 in the center 
to -1-0.25 at 0.5He. 
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ically by +0.20 dex. While this scenario is not plausible 
given the clearly observed color gradients, the test shows 
that the value of the gradient is important. We also con¬ 
sidered a model in which the VTv prior is centered on 
—0.08, which corresponds the mean gradient generally 
found for large samples of ellipticals as discussed in Sec¬ 
tion [5] This more realistic shift has a much milder effect 
on the resulting constraints (see pink star in Figure 14). 

Finally, although ad hoc, it is useful to consider the 
influence of systems where the fit quality is poorer and 
the mass models may be inadequate. The exclusion of 
CSWA6 and CSWA141 has a significant effect on both 
(asps) and (7dm), as described in Section 9.2. 

Systematic uncertainties have a larger effect on asPS 
and 7 dm than on the total density slope 7tot ■ The modi¬ 
fications to our mass models described above shift (ytot) 
by only ±0.07, as quoted in Section 

10. DISCUSSION 

We have leveraged the discovery of large numbers 
of intermediate-separation strong lenses (e.g., Limousin 


et al.||20d9 More et ^|2012 Stark et al.||2013 i to study 

the mass distribution of ETGs central to ~ 10^"* Mq 
group-scale halos. By combining strong lensing, satel¬ 
lite kinematics, and resolved stellar kinematics within 
the lens galaxy, we have constrained the mass distribu¬ 
tion on multiple scales. This is the first study to apply 
these observational methods to a sample of group-scale 
lenses rather than to an individual system. It bridges a 
gap between earlier well-studied samples of lenses at the 
galaxy and cluster scales. Here we consider the impli¬ 
cations of our findings for ETG formation scenarios and 
the dominant physical processes by which DM is redis¬ 
tributed during galaxy formation. 

10 .1. Halo Concentrations at the Group Scale 

On large scales, we have measured the structure of the 
group-scale DM halos via the MGR (Section]^. We find 
a mean concentration that is consistent with expectations 
for unmodified CDM halos, with an uncertainty compa¬ 
rable to the variation among theoretical MCRs based on 
different cosmologies or halo selection criteria. This im¬ 
plies that baryons have minimally affected the structure 
of their halos on these scales, which is in good agreement 
with theoretical expectations but disagrees with some 
earlier observations at this mass scale. On the simula¬ 


tion front, Duffy et al. (|2010| found that the inclusion 
of baryons m group-scale halos increases their concen¬ 
tration by only Ac = 0.1 in their weak feedback sce¬ 
nario, while for stronger feedback prescriptions the con- 
centr ations were low e r than DM-only simulations. Like¬ 
wise, Schaller et al.| (2015a) detect an increase of only 
^C2oo = 0.4 (i.e., 9% or AlogC2oo — 0-04 dex) when 
baryons are i nclu ded in the EAGLE simulations. 

In Section |7.1| we compared our findings with ear¬ 
lier obse rvational constrain ts. Our results are consistent 
with thelBuote et al.| (|2007|) X-ray-based study and with 
Auger et ai.|(|2013|), who studied a lens sample that over- 
laps ours but used different observables. Other lensing 
studies, however, have found higher concentrations at 
halo masses near lO^'^M©. These discrepancies proba¬ 
bly arise from different distributions of 0Ein in samples 
selected in different ways. Samples consisting only of 


lenses with 0Ein > 5" (e.g., Oguri et al. 2012) will nec¬ 
essarily show high concentrations for ~ 10 ^^ Mq halos. 
Lenses with 0Ein — 275 — 5" are much more numer¬ 
ous. Our study focuses on these lenses, which may well 
be the dominant lens population in ~ 10 ^"* Mq halos. 
Tremendous progress has been made by several groups 
(see Introduction) in building up large samples of these 
intermediate-separation lenses. Further progress in using 
them to constrain the MCR will likely require a more 
precise quant ification of the sele ction functions of the 
surveys (e.g., Gavazzi et al.||2014). 


10 .2. Trends in the Total Density Profile Slope and 
Implications for the Bulge-Halo Conspiracy 

A well-known and striking result from lensing and dy¬ 
namical stndies of massive ETGs is the similarity of the 
total density profile to the isothermal form, ptot cc 
which corresponds to a flat circular velocity cur ve (e.g., 
Kronawitter et ar]|2000[ [Gerhard et al.|[2001[ [Treu fcj 


Koopmans 2004|). The proximity of the total density 
slope to the isothermal profile and its relatively small 


halo conspiracy” (e.g., 

Koopmans et al. 

2006 

2009 

Gavazzi et al. 

2UU7I, in which some interaction be 

tween 


the stellar and DiVl components is hypothesized to drive 
the total mass density toward this form, e.g., perhaps as 
fundamental dynamical co nsequence of violent relaxation 
( Treu fc Koopimans||2002[ ) . 

T'hese measurements pertain to ETGs in a fairly nar¬ 
row range of halo masses (ctwmsoo — SLAGS 

galaxy-scale lenses; Section]^. By extending the same 
lensing and dynamical methods to the central galaxies 
of larger halos, we are able to better define the scope of 
the “conspiracy.” First, ETGs are not limited to nearly 
isothermal profiles, but display a wide range of density 


profiles within that spans 7tot — 1.0 — 2.4 (Figure 10 
bottom left panel). Second, 7tot depends jointly on two 
parameters, namely E* and M200, that reflect proper¬ 
ties of the stellar and DM distributions smarately (Fig¬ 
ure 10 bottom right panel, and Equation^. Third, the 
smaTTscatter in 7tot at fixed dEin can be reproduced by 
a range of halo occupation models in which independent 
scatter is assumed in the stellar mass-halo mass, stel¬ 
lar mas s-r adius, and halo mass-concentration relations 
(Figure [T^. Our models without halo contraction explic¬ 
itly have no “cross-talk” between the baryons and DM 
halo, other than what is contained in the stellar mass- 
halo mass relation, yet are able to reproduce the observed 
small scatter in 7tot- Th is third point was also found in 
recent similar models by Dutton & Treu (2014). 

As pointed out by those authors, current cosmologi¬ 
cal, hydrodynamical simulations cannot match both the 
densit y slope and t he star formation efficiency of ETGs 

S et al. 201()[ [Joha nsson et al. 20121 [Dubois et al. 
_I ^ ^museT^ [2018 ), which presumably reflects the 


inadequacy of current feedback implementations. Our 
measurement of the dependence of 7tot on stellar struc¬ 
ture and halo mass (Equation [^ can be used to test the 
performance of different feedback prescriptions in halos 
of widely varying mass. For instance, based on their suite 
of cosmological zoom-in simulations, Remus et al. (2013) 
conclude the isothermal density profile is a “natural at- 
tractor” for the evolution of ETGs and that small devia¬ 
tions from it are compatible with the variation in assem- 
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Fig. 15.— Constraints on 7dm? the average DM density slope 
within Rei for galaxy, group, and cluster-scale ETC lenses. For 
group and cluster lenses, 7 dm is computed following Equation 
where p and M now refer to the DM profile; outer error bars include 
our systematic error estimate added in quadrature. We emphasize 
that 7tot refers to the average slope within Re and cannot be di¬ 
rectly compared to the asymptotic inner slope /d that appears in the 
gNFW profile (Equation!^ and is quoted in other works (including 
N13b). Results f or galaxy-scal e lenses are compiled as indicated in 
the caption. The|Grillo||j2012|| constraint assumes a Salpeter IMF; 
the other results do no t assume a particular IMF a priori. The 
|Sonnenfeld et al.| (|2015|l result is converted from the gNFW inner 
slope to our dennition of 7 d m- L ines follow the halo occupation 
models introduced in Section l9.ll 

bly histories (e.g., number and gas content of mergers) 
seen in the simulations. Their density profiles, however, 
appear to be too homogeneous: nearly all of the simu¬ 
lated ETGs have ytot ^ 1-9, including systems as massive 
as BCGs, which is unlike the broad and systematic trends 
discussed in this paper. A more detailed comparison will 
require a careful matching of the definition of ytot used 
in simulations and observations. 


10.3. The Dark Matter Density Profile within the 
Effective Radius: Galaxies, Groups, and Glusters 

Although the slope of the total density profile clearly 
varies systematically with halo mass, assessing the uni¬ 
versality of the d ark matter profile on small scales is more 
difficult. Figure 15 summarizes the current constraints 
on the inner DM density profile of massive ETGs derived 
from lensing and dynamical observations. 

The situation is most uncertain in the galaxy-scale 
lenses, where the high degree of baryon dominance makes 
it especially difficult to extract detailed information 
about the DM distribution on small scale s (th e “bulge- 
halo degeneracy”). Our analysis in Section 9.1 disfavored 
contra cted models. Sim i lar ana lyses by [Dutton fc Treu| 
(2014) and Auger et al. (2010b I likewise support nearly 
unmodified INE W halos, with mildly contracted ones per¬ 
mitted in the latter study. However, other stud ies in¬ 
stead favor contracted halos at this mass scale. lOgi: 


et al. (2014) incorporate microlensing constraints to in- 
fer a steep DM slope in galaxy-scale lenses, alt hough the 


precision of t he result has been questioned (Schechter 
et al. 2014). Grille (2012) construct a composite den- 


sity prohle based on an ensemble of galaxy-scale lenses 
and, under the assumption of a Salpeter IMF, derive a 


steep slope 7 dm = l-77^og2. The uncertainty, however, 
is large, and Dutton & Treu (2014[) revised this analysis 
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Sonnenfeld et al. 


( 2012 ) 


to a shallower slope of 1.40^J 
analyzed a rare double Einstein' ring, allowing them to 
break the usual bulge-halo degeneracy and measure a 
slope 7 dm = 1-7 ± 0.2 that is indicative of a contracted 
halo, but it is unclear how representative this single lens 
is of the overall p opulation of massive ETGs. Building 


on the meth ods of Treu & Koopmans (2004), Sonnenfeld 
et al. (2014) studied an ensemble of galaxy-scale lenses 


in the SLAOS and SL2S surveys; although the DM slope 
cannot be constrained for individual objects with unre¬ 
solved stellar kinematics, variations in the source red- 
shift lead to variations in i?Ein/7?e, which allows loose 
constraints to be placed on the DM slope. Their results 
favor NEW or even shallower DM profiles, although the 
authors stress the large uncertainties. 

Greater progress has been made in constraining the 
DM profile in higher-mass halos owing to the reduced 
dominance of baryons in these systems and the richer 
observational constraints. N13a,b showed that the cen¬ 
tral galaxies in massive clusters have a DM profile that is 
shallower than the canonical NEW slope within ~ i?e, a 
conclusion that is further supported by the halo occupa¬ 
tion models presented here. Several observations guided 
the N13 interpretation of this result: the similarity of ra¬ 
dial scale at which the DM profile deviates to the effective 
radius of the BGG, the remarkable similarity of the total 
density profile to the NEW slope on these scales, and a 
tentative correlation between the inner DM profile and 
i?e when individual objects are examined. These lines 
of evidence suggest that the flattening the DM profile 
within BGGs is connected with the assembly of stars in 
the galaxy. Such a connection might be achieved dynam¬ 
ically through the transfer of angular momentum from 
infalling satellites to the DM halo via dynamical fric- 


tion ( 

El-Zant et al.|20011 2004 

Nipoti et al.|2004| Tonini 

et al. 

2U(J6l IKomano-Diaz et a 

1.1120081 jJardel & Sellwood 


Johansson et al. 2009 Del Fopolo 2012). Laporte 

& White (2U15|) recently used iV-body simulations con- 


mix during the assembly of a BCG. They verified that 
the DM profile can be flattened over time, provided that 
the assembly of the BCG is dominated by dissipationless 
stellar accretion following 2 ~ 2 and that the progenitor 
galaxies at z ~ 2 are realistic. 

Our new observations of group-scale lenses, located in 
halos that are ~ 10 times less massive than the N13 clus¬ 
ters, instead favor unmod ified NEW halos. M ildly con¬ 
tracted halos following the Gnedin et al. (2004) model are 
also allowed at the 1.4 —2 .Oct level, but stronger contrac¬ 
tion prescriptions are ruled out if the stellar mass scale 
is indeed close to that of a Salpeter IMF. Cored DM ha¬ 
los with rcore/?'200 as large as those seen in clusters are 
also disfavored by our data. They are not yet definitively 
ruled out: Figure [T^ shows that the groups deviate from 
the cored halo model by < 2cr if CSWA6 and CSWA141 
(whose velocity dispersion profiles are less well fit) are 
excluded from the present sample. 


10.4. Interpreting Possible Trends in the Inner Dark 
Matter Density Profile 
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These observations may support a scenario in which 
the DM profile is determined by the relative impor tance 


of “two pha ses” of massive galaxy formation (e.g., Oser 


et al. 2010): early compression of the halo is associated 


with dissipative, in situ star formation fed by gas infall, 
while later reduction of the central density through dy¬ 
namical friction occurs as dissipationless merging domi¬ 
nates the mass assembly. The balance of these processes, 
given by the fraction of stars formed in situ , is a strong 
function of halo mass: Behroozi et al. (2013) find roughly 
46%, 24%, and 7% at z = 0.3 for halos that correspond 
to our galaxy, group, and cluster lens samples. This may 
well lead to systematic trends in the inner DM profile as a 
function of halo mass, with the underlying physical driver 
being the variation in star formation efficiency M*/M 2 oo 
and assembly history, or in other words, the balance of 


dissipational a nd dissipationless formation (Lackner & 
Ostriker||2010 1. 


In this picture halo contraction becomes increasingly 
significant for ETGs in lower-mass halos, where more 
gas has cooled and been transformed into stars, whereas 
flattened DM profiles appear at the highest masses due 
to the rising importance of dissipationless mergers. Such 
a trend could explain why no single model in Figure 
is able to fit observations from galaxy to cluster scales. 
It is also possible that AGN feedback plays an increasing 
role in flattening the DM cusp in more massive systems 
( Martizzi et al.|2013 ). This paper is an early step toward 
quantifying these trends observationally. 

Although it might be supposed that absence of shal¬ 
low DM profiles in at least some halos is a blow to in¬ 
terpretations that invoke the particle microphysics, such 


as se lf-interactions (e.g., Peter et al. 2013 


Rocha et al. 


2013), this is not necessarily the cas^ Appealing to a 


velocity-dependent scattering cross-section ctdm does not 
help if udm and v vary inversely, as is usually presumed: 
we find stronger evidence for cores in higher mass halos, 
which would have lower (Tdm and smaller cores. Instead, 
the main uncertainty is the behavior of self-interacting 
dark matter in a realistic setting with baryons. Further 
simulations are needed to test this. 


Schaller et al. (2015b) recently compared the N13 clus- 
ter observations with their FAGLF simulations. Unlike 


the simulations by Martizzi et al. (2013) and Laporte 


& White (2015), they found no evidence tor shaifow DM 


sfopes and suggested instead that neglecting radial trends 
in velocity anisotropy could bias the observations. The 
volume of their cosmological simulation was too small 
to contain very massive systems analogous to the N13 
clusters, so they selected 6 halos with masses intermedi¬ 
ate between the N13 sample and the present group-scale 
lenses. Since we do find evidence for cores in 10^"* M© 
halos, we suggest that a mismatch of hal o masses could 
be partially responsible. In Section |9.3| we considered 
a model with radially varying anisotropy and found the 
results to be very comparable to the range of constant 
anisotropy models that were considered by N13. This 
is reassuring, although a wider range for the anisotropy 
profile could be considered. Further progress can be 
made with future observations using the new integral 
field spectrographs appearing on large telescopes (e.g., 
KCWI and MUSF), which will allow higher-order mo¬ 
ments of the velocity distribution to be measured at large 
radii in these lenses. 


10.5. Future Improvements 

Larger samples of lenses with data of similar quality 
are needed to confirm the possible difference we find be¬ 
tween the inner DM profiles of groups and clusters. Fur¬ 
ther progress will also come from improved understand¬ 
ing of the stellar populations. Some of our results for the 
group-scale lenses are sensitive to the radial gradient in 
M*/L. This probably plays a lesser role in galaxy-scale 
lenses, which typically do not have resolved stellar kine¬ 
matic data, and also in clusters, where we have verified 
that the BGGs in the N13 sample show weak color gra¬ 
dients (consistent with a mean of zero) over the relevant 
radial range. However, it has a surprisingly strong effect 
in the BGGs. Age and metallicity gradients could be bet¬ 
ter tested with near-infrared photometry, but it is also 
possible that the IMF itself varies with radius. Testing 
this requires high-quality spectroscopy at large radii, and 
currently we have littl e data to constrain the magn itude 
of such a gradient (e.g., Martm-Navarro et al. 2015). For¬ 
tunately, improved constraints are within reach of future 
observational programs. 


11. SUMMARY 

1. We assembled a sample of 10 group-scale lenses 
(^^Ein = 2.5"—5.1") with high-quality strong lensing 
and kinematic data, including new Keck observa¬ 
tions of 8 systems. This is the largest uniformly 
analyzed sample constructed with the aim of mea¬ 
suring the mass distribution over a wide range of 
scales, from within the central galaxy to the virial 
radius. 

2. Measurements of the velocity dispersion of satellite 
galaxies indicate a mean halo mass of (log M 200 ) = 
14.0 ± 0.1, intermediate between earlier samples of 
galaxy- and cluster-scale lenses. 

3. The mean concentration of the DM halos is 
(c 20 o) = 5.0 ± 0.8. Correcting for the expected 
bias arising from our selection in Einstein radius, 
relative to the underlying ETG population, yields 
(c 2 oo) = 3.9 ± 0.6 (Section]^. 

4. These concentrations are consistent with the cur¬ 
rent range of theoretical MCRs, suggesting that 
baryons have a minor effect on the structure of their 
halos at large radii near the scale radius. This is 
consistent with hydrodynamical simula tions (e.g.. 


Duffy et al. 

2010 

Schaller et al.||2015a|); some ear- 

lier observa 
was likely a 

dona) 

work, e.g., 

Uguri et al. (|2012), 

ttected by incomplete treatment ot how 


strong lenses were selected. 


5. The average slope of the total density profile within 
Re is (7tot) = 1.64±0.05 (stat.)±0.07 (sys.) for our 
group-scale lenses (Section |^. By combining with 
earlier lens samples to examine trends over a factor 
of ~ 60 in halo mass, we show that ytot depends 
on both the stellar surface density E*, as shown in 
earlier lensing studies, but also on the halo mass. 
This reflects the strong non-homology of massive 
ETGs and shows that their internal density struc¬ 
ture cannot be derived from the stellar light alone. 
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6. We built a set of simple halo occupation models 
based on the M* — M 200 and M* — Re relat ions ob¬ 
served in our lens compilation (Section 9.1). These 
reproduce earlier findings that galaxy-scale lenses 
favor a Salpeter-type normalization of the IMF 
(e.g., Auger et al.|2010b Treu et al.|2010 ) and that 
the DM prohle at the centers of m assive clusters is 


shallowe r than the NFW profile (|Newman et al. 
2013a|bD . 


7. Observations of the group lenses alone do not ro¬ 
bustly resolve the covariance between M*/L (and 
therefore the IMF) and the inner DM profile. How¬ 
ever, if the Salpeter-type IMF inferred for the 
galaxy- and cluster-scale lenses holds for the cen¬ 
tral galaxies of groups—as might be expected from 
the similarity of their central velocity dispersions— 
then the group lenses favor unmodified or weakly 
contracted halos. They disfavor an inner DM den¬ 
sity profile as shallow as those found for clusters, 
although this difference is not yet definitive. 


8. Combining our results with earlier studies, we sug¬ 
gest that the effect of baryons on the structure of 
their DM halos may vary with halo mass due to un¬ 


derlying trends with star formation efficiency and 
assembly history. Clusters show evidence for a re¬ 
duction in the central DM density, which may arise 
from the dominant role of dissipationless merging 
in the assembly of BCGs. In group- and galaxy- 
scale halos, more of the gas has cooled and been 
transformed into stars, leading to a greater role for 
dissipation and adiabatic contraction. Larger lens 
samples with high-quality data are required to test 
this tentative variation in the small-scale distribu¬ 
tion of DM. 
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APPENDIX 

A. NOTES ON INDIVIDUAL LENS MODELS 

CSWA107 has a naked-cusp configuration. Not much structure is visible in the current ground-based imaging, so 
we do not include external shear in this model. As for the other three-image lenses, we fix the center of mass to the 
light but allow its ellipticity and position angle to vary independently. 

CSWA141 is the most ambiguous system. The bright blue clump is clearly singly imaged, while its fainter, extended 
tail may be composed of merging multiple images of a fainter part of the source. Furthermore, a second galaxy 
contributes significantly to the lensing potential. Given the increased uncertainties in this system, we hx the center, 
ellipticity, and position angle of the main deflector to follow the BCG light and omit external shear. We also increase 
the uncertainty in k to 0.1. 

CSWA164 is a nearly complete Ei nstein ring. Our measuremen t of 6>q in = 3768 and the low external shear F = 0.022 
agree closely with measurements by Kostrzewa-Rutkowska et al. (20141. 

CSWA165 has a qua d conhguration that constrains the lens model well. Our 0Ein = 4733 is slightly larger than the 


Kostrzewa-Rutkowska et al. 


3777to;i6 measured by 
depth and resolution than our DEIMOS i maging 


(2014), but their analysis was based on SDSS images with poorer 


CSWA6 (The Clone) was di scovered by L in et a l. (2009), who modeled the deflector using a single SIE with 0Ein = 
3782, smaller than our 4736. Jones et al. 1 2010 FTAowever, analyzed HST imaging and followed a procedure more 
analogous to our own, separating the contribution of the satellite galaxies. They measure Einstein radii of 370 ± 0'.'3, 
0'.'5±073, and 0'.'3±072 for their Cl, G2, and G3 deflectors, where G1 represents the BGG and G2 and G3 are satellites. 
These agree well with our measurements of 372, 074, and 075, respectively. Satellites in our model contribute 21% of 
the mass within 0Ein, larger than for the other lenses. Owing to the increased uncertainty in decomposing the mass 
components, we increase the error in k to 0.1. 

GSWA7 forms multiple images of two background galaxies, each presenting 3 images visible in the HST data. Given 
the presence of the multiple sources and the compactness of their images, we decided to use glafic (|Oguri||2010| 
to model this lens using the positions of the six images as constraints. This accounts for the different appearance of 
GSWA7 in FigureAs in the other naked cusp lenses, we fix the center of mass to that of the BGG. We hnd a fairly 
high external shear of F = 0.149 consistent with the prominent background structure seen in our redshift survey. The 
axis ratios of the mass and light agr e e (gma ss = 0.67, gught = 0.65), suggesting that the ellipticity and shear have been 
successfully separated. Kubo et al. (2009) noted that the arc spectrum shows superposed spectra at z = 1.411 and 
z = 1.38, and we see the same m our DEIMOS spectrum. By fixing the redshift of the outer brighter system to 1.411 
and allowing that of the fainter inner system to vary in our lens model, we find a redshift of 1.38. This confirms that 
the fainter multiply imaged system is the galaxy seen in absorption in the spectrum of the brighter background galaxy. 
Our 0Ein = 2773, which refers to zs = 1.411, is smaller than the 377 estimated in the discovery paper by Kubo et ah, 
but this can likely be attributed to their neglec t of the (high) ellip ticity. 

EOGL (Eight O’Glock Arc) was discovered by|Allam et al. (|20 07t, who reported a preliminary S IE-base d 6>Ein = 2791. 
Upda ted analyses with HST observations revised this to 3. 32 (|Dessauges-Zavadsky et al.||2011l see also|Shirazi et al. 
20T4|, very close to our 3729. 

JI59413-1100 has the most regular light distribution of the 13 group-scale lenses in the SL2S sample presented by 
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Limousin et al. (2009), suggestive of a relaxed system. We find that external shear is essential to accurately fit the 
features seen in the HST imaging. Reassuringly, when external shear is incorporated into the model, the axis ratio of 
the mass distribution (g^ass = 0.55) closely matches that of the BGG light at the arc radius. (The overall axis ratio 
of the BGG in Table [mis higher due to the presence of an ellipticity gradient.) This arc is the only one in our sample 
that lacks a spectroscopic redshift. We detect its continuum throughout the DEIMOS wavelength range, but despite 
its blue color see no emission lines. Based on the non-detect ion of [O II), we plac e a lower limit z > 1.6. Using the 
uriz colors measured in public CFHT imaging, we use EAZY (Brammer et al.||2008) to infer a 99% upper limit on the 
photometric redshift of z < 2.45. This reflects the lack of a break in the u photometry that would occur at higher 
redshifts due to absorption by the intergalactic medium. We use a fiducial zs = 2.0 in our analysis and note that the 
leasing distance ratio D^s/Dg changes by only ^ ±5% over the allowed range zs = 1.6 — 2.45. We consider this as an 
additional random error in R that is added in quadrature. 

CSWA163 is constrained only by SDSS imaging, w hich has relatively po or quality. We therefore omit external shear 
from this model. Our dEin = 3"49 agrees closely with||Deason et aT (2013). 


CSWAl ( The Cosmic Horseshoe) has the largest Einstein radius in our sample and was studied in detail by Dye 


et al. 
we fini 


2008|). The 0Ein = 4"97 obtained for their SIE model is within 2% of our 5708. In agreement with Dye et al. 
a very low contribution from external shear (T = 0.022). 


B. HIERARCHICAL BAYESIAN INFERENCE METHOD 

Here we describe the mathematical framework that we used to analyze the MGR. In this approach, we assume that 
the group lenses are drawn from a parent sample of lenses which is characterized by independent Gaussian distributions 
of logM 2 oo and log c%n = logc^nn + 0.07(log M^nn — 14). Here we account for the slope of MGR using the theoretically 
expected slope (Duffy et al.||2008), since the range of masses in our sample is much too narrow to constrain it. With 
the slope of the MGR removed, we may consider \ogM 200 and logc^oo to be independently distributed. 

Our aim is to constrain the mean and intrinsic dispersion of this parent distribution—(logM 2 oo)j (logc^oo)) criogM, 
and tJiog c') which are known as the hyperparameters—using the posterior distributions of each of the 8 lenses (those 
with M200 measured from satellite dynamics) in the (M200, C 200 ) plane, after marginalizing over the other parameters. 
For brevity, we denote the hyperparameters as (M), (c), ctm, and (7^, respectively, and the set of these as w. Their 
posterior distribution is then 


P{uj\D) oc Pr{uj)P{D\u}) 

= Pr(a;) JJ dMdc Pi{Di\M, c)N{M, c;u}) 

ocPr(w)n // dM (M, c; uj). 


(Bl) 


Here Di represents the data for lens i and D the combined data set. N{oj) is the trial parent distribution, i.e., a 
bivariate Gaussian with mean ((M), (c)) and dispersion (ctm, cTc)- The posterior distribution Pi{M, c\Di) for each lens 
is constructed from its MCMC chains. Since both the prior Pri{M, c) used i n th e first-level inference for each lens and 


the prior Pr(uj) on the hyperparameters lo are broad and uniform. Equation Bl amounts to a product of i integrals of 
the trial parent distribution m ultiplied by the posterior for lens i. We sample P{lo\D) using the MCMC code emcee 
( Foreman-Mackey et al.|[2013 ). 


Similar techniques have been used by, e.g. 
apply the same framework to study the distributions of 7 dm and asPS ■ 


Bolton et al. 

(2012), N13a, and Auger et al. 

(2013 

). In Section 9.2 


we 
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